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Summary 

This  investigation  concerns  equilibrium  fields  with  discontinuous 
displacement  gradients,  but  continuous  displacements,  in  the  theory  of  fi- 
nite plane  deformations  of  possibly  anisotropic,  compressible  elastic  sol- 
ids. "Elastostatic  shocks"  of  this  kind,  which  resemble  in  many  respects 
gas-dynamical  shocks  associated  with  steady  flows,  are  shown  to  exist  on- 
ly if  and  when  the  governing  field  equations  of  equilibrium  suffer  a loss  of 
ellipticity.  The  local  structure  of  such  shocks,  near  a point  on  the  shock- 
line, is  studied  with  particular  attention  to  weak  shocks,  and  an  example 
pertaining  to  a sho’dk  of  finite  strength  is  explored  in  detail.  Also,  neces- 
sary and  sufficient  conditions  for  the  "dissipativity"  of  time -dependent 
equilibrium  shocks  are  established.  Finally,  the  relevance  of  the  analysis 
carried  out  here  to  localized  shear  failures  — such  as  those  involved  in  the 
formation  of  Liiders  bands  — is  discussed. 


Introduction 


Several  years  ago  — in  connection  with  asymptotic  studies  of 


*The  results  communicated  in  this  paper  were  obtained  in  the  course  of  an 
investigation  supported  by  Contract  N00014-75-C-0196  with  the  Office  of 
Naval  Research  in  Washington,  D.  C. 
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crack  problems  within  finite  elasticity  theory  — we  encountered  the  fact 
that  the  field  equations  of  nonlinear  elastostatics  may  suffer  a loss  of  ellip- 
ticity  for  some  elastic  materials  in  the  presence  of  sufficiently  severe  lo- 
cal deformations.  A special  isotropic  material  of  this  kind  was  explored 
in  [l],  where  we  determined  all  three-dimensional  homogeneous  deforma- 
tions at  which  such  a breakdown  of  ellipticity  occurs.  More  recently  [z], 
we  deduced  explicit  necessary  and  sufficient  conditions,  in  terms  of  re- 
strictions upon  the  principal  stretches  and  upon  the  dependence  of  the 
strain-energy  density  on  the  principal  stretches,  for  ordinary  and  strong 
ellipticity  of  the  equations  governing  finite  plane  deformations  of  general 
homogeneous  and  isotropic,  compressible  hyperelastic  solids. 

A failure  of  ellipticity  of  the  elastostatic  equations  appropriate  to 
such  solids  suggests  the  possible  emergence  of  solution  fields  exhibiting  a 
loss  of  smoothness.  In  the  present  paper  our  interest  is  directed  in  par- 
ticular at  the  possibility  of  solutions  that  possess  finite  jump  discontinu- 
ities in  the  first  displacement  gradients  across  certain  curves,  while  the 
displacements  themselves  still  remain  everywhere  continuous.  We  call 
solutions  of  this  kind  "elastostatic  shocks"  in  allusion  to  gas-dynamical 
shocks  associated  with  stationary  inviscid  flows,  to  which  they  bear  a 
more  than  casual  resemblance.  Moreover,  since  our  current  concern  is 
with  the  local  state  of  affairs  near  an  interior  point  on  a shock-line  (line  of 
displacement-gradient  discontinuity),  it  will  be  sufficient  to  confine  our  at- 
tention here  to  piecewise  homogeneous  elastostatic  shocks:  thus  we  assume 
the  shock-line  to  be  a straight  line,  on  either  side  of  which  the  plane  defor- 
mation under  consideration  is  homogeneous. 

Before  we  proceed  to  outline  the  content  of  this  paper,  we  should 
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make  clear  that  our  motivation  in  studying  elastostatic  shocks  is  physical 
in  origin.  Discontinuities  of  the  kind  described  above  arise  typically  in 
the  idealization  of  localized  shear  failures,  such  as  those  involved  in  the 
formation  of  Liiders  bands. 

In  Section  1 we  recall  some  prerequisites  from  the  nonlinear  theory 
of  plane  elastostatic  deformations,  including  the  relevant  notions  of  ellip- 
ticity.  Here  we  also  cite  the  explicit  conditions  of  ordinary  and  strong  el- 
lipticity  established  in  [2]  for  the  special  case  of  material  isotropy,  as 
well  as  results  derived  there  regarding  the  inclination  of  the  characteris- 
tic curves  that  accompany  a loss  of  ordinary  ellipticity. 

Section  2 contains  a detailed  discussion  of  response  properties  in 
plane  strain  of  two  special  isotropic  elastic  materials,  which  are  used  in 
the  remainder  of  the  paper  to  illustrate  various  general  conclusions.  The 
first  of  these  materials,  which  was  also  the  object  of  [ 1 ] , is  governed  by 
a stored-energy  function  originally  proposed  by  Blatz  and  Ko  [3]  on  the  ba- 
sis of  experiments  with  a highly  deformable  rubberlike  material.  The  sec- 
ond particular  material  introduced  in  Section  2 is  strictly  hypothetical  iri 
nature  and  is  employed  later  on  to  bring  into  evidence  certain  constitution- 
ally determined  qualitative  differences  in  behavior. 

A complete  definition  of  piecewise  homogeneous  elastostatic  shocks 
is  spelled  out  in  Section  3,  where  their  kinematics  is  analyzed  in  detail. 
Further,  we  prove  here  that  the  existence  of  such  an  equilibrium  shock  is 
contingent  upon  a loss  of  strong  ellipticity  of  the  displacement  equations  of 
equilibrium,  appropriate  to  the  material  at  hand,  at  some  homogeneous  de- 
formation. The  results  obtained  in  Section  3 — like  most  subsequent  re- 
sults — are  not  restricted  to  isotropic  materials. 


r 
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The  considerations  carried  out  in  Section  3 suggest  that  if  there  ex- 
ists any  piecewise  homogeneous  shock  corresponding  to  a pre-assigned, 
deformation  on  one  side  of  the  shock-line,  there  is  a one-parameter  fami- 
ly of  such  shocks.  In  Section  4 we  take  for  granted  the  existence  of  a 
smooth  family  of  shocks  in  a neighborhood  of  the  shockless  state.  We  then 
deduce  various  "weak-shock”  results  — encompassing  the  jumps  across 
the  shock-line  of  physically  significant  field  quantities  — to  dominant  or- 
der in  a shock-strength  parameter  that  measures  the  departure  from  the 
shockless  solution.  One  of  the  conclusions  reached  in  this  manner  is  that 
the  existence  of  the  presupposed  family  of  equilibrium  shocks  demands  a 
loss  of  ordinary  ellipticity  at  the  given  homogeneous  deformation.  More- 
over, the  limiting  shock-line  at  zero  shock- strength  is  found  to  be  a char- 
acteristic line  associated  with  this  deformation.  It  should  be  mentioned 
that  Rudnicki  and  Rice  [4]^  had  previously  arrived  at  the  required  failure 
of  ellipticity  in  dealing  with  weak  shocks  of  a related  type  for  a broader 
class  of  materials  that  includes  the  elastic  solid. 

Section  5 is  devoted  to  an  instructive  example  of  a global  analysis 
of  a piecewise  homogeneous  equilibrium  shock  of  finite  strength,  based  on  • 

the  first  of  the  two  special  isotropic  elastic  materials  discussed  in  Sec- 
tion 2.  In  this  instance  explicit  results  for  all  possible  shocks  can  be  ob- 
tained in  a transparent  elementary  form.  In  particular  we  find  here  that 
every  shock  admitted  is  accompanied  by  a breakdown  of  ordinary  elliptic- 
ity at  least  on  one  side  of  the  shock-line,  regardless  of  the  strength  of  the 
shock. 

In  gas  dynamics  the  sign  f *■  the  shock  strength  is  determined  by  the 
^See  also  Rice's  [ 5]  more  recent  paper. 


entropy  inequality,  which  models  the  dissipative  character  of  the  process 
of  shock  formation  in  the  absence  of  viscosity.  Guided  by  this  fact  we 
generalize  in  Section  6 the  results  of  Section  3 to  time -dependent  piece- 
wise  homogeneous  quasi- static  shocks  and  then  establish  an  energy  identi- 
ty that  serves  as  a basis  for  a proposed  criterion  of  dissipativity  appropri- 
ate to  such  shocks.  Thereafter  we  obtain  necessary  and  sufficient  dissipa- 
tivity conditions,  which  lead  to  an  inequality  analogous  to  the  entropy  con- 
dition of  gas-dynamics.  The  dissipation  inequality  is  applied  at  the  end  of 
Section  6 to  determine  the  sign  of  the  shock  strength  in  the  special  global 
solution  arrived  at  in  Section  5,  as  well  as  for  weak  shocks  in  the  second 
of  the  two  particular  isotropic  materials  mentioned  earlier.  This  sign,  in 
turn,  governs  the  sign  of  the  jump  in  the  mass  density  as  the  shock-line  is 
traversed.  We  find  that  the  shock  strength  may  be  positive  or  negative  de- 
pending On  the  particular  nature  of  the  elastic  material  experiencing  the 
shock. 

In  Section  7 we  combine  the  special  results  of  Section  5 and  the  gen- 
eral conclusions  reached  in  Section  6 in  ein  attempt  to  illustrate  the  rele- 
vance of  elastostatic  shocks  to  the  phenomenon  of  Ltiders  bands  in  a slab 
under  uni-axial  tension  or  compression.  For  this  purpose  we  view  the 
evolution  of  such  a band  — in  the  vicinity  of  a point  on  the  interface  be- 
tween band  and  slab  — as  a bifurcation  from  a homogeneous  deformation 
that  has  become  dynamically  unstable  following  a loss  of  ellipticity,  into  a 
time -dependent  piecewise  homogeneous  equilibrium  shock.  The  ensuing 
results,  though  tied  to  a very  limited  elastic  material,  display  various 
striking  qualitative  features  that  are  reminiscent  of  familiar  experimental 
observations.  In  this  connection  we  refer  to  the  work  of  Hill  and 


Hutchinson  [6],  where  various  bifurcations  from  a state  of  plane-strain  uni- 
axial tension  are  investigated  for  a class  of  incompressible,  incrementally 
linear  time- independent  materials.  Among  these  bifurcations  are  localized 
shearing  modes. 

1 , Preliminaries  on  finite  elastostatic  plane  strain. 

In  this  expository  section  we  recall  from  [2]^  certain  results  pertain- 
ing to  the  nonlinear  theory  of  plane  elastostatic  deformations  which  are  es- 
sential to  the  analysis  of  the  class  of  discontinuous  plane  elastostatic  fields 
that  constitutes  our  main  objective. 

As  for  notation,  we  shall  use  boldface  letters  to  denote  vectors  and 
second-order  tensors  in  two  dimensions,  as  well  as  two-rowed  column  and 
square  matrices.  Further,  the  same  boldface  letter  will  be  employed  to 
designate  a vector  or  tensor  and  its  matrix  of  scalar  components  in  the  un- 
derlying rectangular  cartesian  coordinate  frame. 

Let  ft  be  the  open  region  of  the  (x^ , x^) -plane  occupied  by  the  interior 
of  the  middle  cross-section  of  a cylindrical  or  prismatic  body  in  its  unde- 
formed configuration.  A plane  deformation  of  such  a body  — parallel  to  the 
(xpX^) -plane  — is  described  by  a suitably  smooth  and  invertible  transfor- 
mation 


j^=^(x)  = x + u(x)  for  all  x in  ft. 


(1.  1) 


which  maps  ft  onto  a domain  ft^  of  the  same  plane.  Here  x is  the  position 
vector  of  a generic  point  in  R,^(x)  is  its  deformation  image  in  ft^^.,  while  u 
is  the  displacement  vector  field.  Accordingly,  x^  and  y^  are  the  cartesian 

^ Parts  of  this  resume^  are  taken  verbatim  from  Section  1 of  [2]. 
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rnaterial  and  spatial  coordinates,  respectively.  We  call  F the  deformation- 
gradient  tensor  belonging  to  the  mapping  (1.  1)  and  J its  Jacobian  determi- 
nant (area  ratio).  Thus, 


]f^=V^=[ay^/axp]  , J=detF>Oonft. 


(1.  2) 


Let  T be  the  in-plane  Cauchy  stress -tensor  field  accompanying  the 
deformation  at  hand,  so  that  stands  for  the  prevailing  components  of 
"actual"  or  "true"  stress.  The  appropriate  two-dimensional  stress  equa- 
tions of  equilibrium,  in  the  absence  of  body  forces,  then  take  the  form 


divT=0  , T=t'^  or  9T^p/ayp-0  , on  R,.. 


(1.  3)- 


Next,  suppose  a represents  the  in-plane  Piola  stress-tensor  field  corres- 


ponding to  JT,  whence 


a=JT(F-^)^or 


(1.  4) 


where  a o components  of  ’’nominal”  or  "pseudo  - stres s and  F , 

Qfp 

with  the  components  F'^,  designates  the  inverse  of  the  nonsingular  tensor 
F.  Equations  (1.  2),  (1.  3),  (1.  4)  lead  to  the  equilibrium  conditions 


diva=0  or  aa^p/9xp=0  on  R, 


(1.  5) 


but  a is  in  general  not  symmetric. 


^Greek  subscripts  have  the  range  (1,  2)  and  we  shall  employ  the  usual  sum- 
mation convention. 

^If  M is  a two-by-two  matrix  with  elements  alternatively  write 

[M^^]  in  place  of  M. 


A superscript  T will  always  indicate  transposition. 


Turning  now  to  the  governing  constitutive  relations,  we  assume  that 
the  body  is  elastic  and  possesses  an  elastic  potential  W,  which  represents 


the  strain-energy  density  per  unit  undeformed  volume.  Moreover,  we  re- 
strict the  material  to  be  homogeneous.  Consequently,  W depends  upon  po- 
sition on  ft  exclusively  through  the  deformation-gradient  tensor  F,  and  the 
constitutive  law  — as  far  as  the  in-plane  Piola  stresses  are  concerned  — 
becomes 

a=W„  or  a ,=aW/aF  (1.6) 

~ F ap  ap 


Substituting  from  (1.6)  into  (1 . 5),  and  invoking  (1.1),  (1.2),  one  arrives  at 
the  displacement  equations  of  equilibrium 


c „ c (F)u  =0  on  ft, 
aPy6'~'  Y . 6P 


(1.  V 


provided 


c . JF)=a  W/9F  ,,aF  .. 

apy6~  Y6  »p 


(1.8) 


For  all  unit-vectors  N,  let  Q(N;F)  be  the  symmetric  tensor  defined  by 

The  quasi-linear  system  of  partial  differential  equations  (1. 7)  is  elliptic  at 
a solution  u (with  continuous  first  and  piecewise  continuous  second  partial 
derivatives  on  ft)  and  at  a material  point  x,  provided 

detQ(N;F(x));tO  for  all  with  1n1=1,  (1.  10) 

where  F(x)  is  the  value  at  x of  the  deformation  gradient  field  gene  rated  by  u. 

Subscripts  preceded  by  a comma  indicate  partial  differentiation  with  re- 
spect to  the  corresponding  material  cartesian  coordinates. 

2 

The  symmetry  of  Q(N;F)  follows  from  ‘^apY6~*^Y6<^P" 
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If  (1.  10)  fails  to  hold,  so  that 

detQ(N;F(x))=0  for  some  N with  |n|=1,  (1.11) 

then  N is  normal  to  a characteristic  curve  (referred  to  the  undeformed 
state)  associated  with  the  solution  u at  x.  The  characteristic  curves,  in 
turn,  are  the  only  possible  carriers  of  "weak  discontinuities"  of  u:  across 
such  a curve  the  second  normal  derivative  of  ju  may  exhibit  a finite  jump 
discontinuity,  whereas  its  second  tangential  derivative,  as  well  as  the 
) mixed  tangential -no rmal  derivative  are  bound  to  remain  continuous.  Thus. 

the  ellipticity  condition  (1.  10)  precludes  the  existence  of  real  characteris- 
tics; (1.  10)  is  necessary  and  sufficient  in  order  that  every  solution  jr  of  the 
presupposed  smoothness  be  free  of  weak  discontinuities  and  hence  in  fact 
‘ twice  continuously  differentiable  at  the  material  point  under  consideration. 

Following  common  usage,  we  call  the  system  (1.7)  strongly  elliptic 
at  a solution  .u  and  a material  point  x,  provided  Q(K>E(~))  positive-defi- 
nite for  every  unit-vector  N,  i.  e.  provided 

M*Q{N;F(x))M>0  (1.  12) 

■ for  all  unit-vectors  M and  N,  Clearly,  the  strong  ellipticity  of  (1.  7)  at  u 

and  X implies  its  ordinary  ellipticity. 

We  proceed  next  to  the  particular  case  of  an  isotropic  body  under- 
^ going  a plane  deformation  of  the  form  (1.  1).  To  this  end  let  C and  g,  re- 

spectively, be  the  right  and  the  left  Cauchy-Green  deformation  tensors  as- 
sociated with  the  mapping  (1.  1),  whence 

\ 

C=F^F  , G=FF^.  (1.  13) 
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These  two  symmetric,  positive-definite  tensors  have  t.ie  same  fundamental 

scalar  invariants  and  hence  common  principal  values,  which  are  the  squares 

of  the  local  principal  stretches  (length-ratios);  the  latter  will  be  denoted  by 

\ >0.  In  view  of  (1.  13)  and  the  second  of  (1.2),  the  deformation  invariants 
a 


just  mentioned  obey 


I=trC=trG=F^pF^^=X^  + \2. 


J = ydetC  =VdetG  = detF=  jE^pE^^F^^F 


(1.  14)' 


in  which  stands  for  the  components  of  the  two-dimensional  alternator. 

In  the  special  instance  of  an  isotropic  material  subjected  to  a plane 
deformation  (1.  1),  the  strain-energy  density  W involves  F merely  through 
the  two  deformation  invariants  I and  J.  In  these  circumstances  one  has 


(1.  15) 


where  — in  order  to  avoid  unduly  cumbersome  notation  — we  have  employed 

the  same  functional  symbol  in  three  distinct  connotations.  In  particular, 

2 

W(I,  J)  and  W(Xj,\2)  stand  for  the  plane- strain  elastic  potential  of  the  mate- 
rial at  hand  in  terms  of  the  deformation  invariants  and  the  principal  stretches, 
respectively.  From  (1.  14)  one  has 


^For  convenience  we  use  J,  rather  than  the  common  determinant  of  C and  G, 
as  the  second  invariant. 

2 

These  two  alternative  forms  of  the  potential  are  the  restrictions  to  plane 
deformations  of  their  spatial  counterparts,  which  in  general  depend  upon  the 
three  invariants  of  the  three-dimensional  deformation  tensors  and  on  the 
triplet  of  principal  stretches.  While  W(I,J)  fully  characterizes  the  in-plane 
response  to  a plane  deformation,  it  is  insufficient  for  the  determination  of 
the  out-of-plane  stresses  so  induced. 


t - 
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91  _■?  p 9 J -p  _jp-l 

9F^p"^-^ap'  9F^p"®Py  “6^6y  Po' 


and  (1.  6),  (1.  15),  (1.  16)  furnish  the  stress -deformation  relation! 


a=W„=2W^F+JWT(F'S'^ 


a =2W  F^.+W^So  e ,F. 

oP  9F  „ 1 QfP  J Py  6y 


(1.  16) 


(1.  17) 


in  case  of  material  isotropy,  where  W^.  and  Wj  are  the  partial  derivatives 
of  W(I,J)  with  respect  to  its  arguments.  The  corresponding  constitutive  re- 
lations for  the  actual  stresses  are  obtained  by  means  of  (1. 4)  and,  because 
of  (1.  13),  may  be  written  as 


t='t'^  = ^W^G+W.1 


'^op-'^pa"  J^I*^ap^^J^ap  ’ Sp'^ay^PY’ 


(1.  18) 


if  1 is  the  two-dimensional  idem-tensor  and  6^^  the  Kronecker-delta.  Ac- 
cording to  (1.  18)  the  tensors  _T  and  G have  the  same  principal  axes.  Finally, 
(1.  17)  together  with 


(1.  19) 


yield  the  appropriate  specialization  of  (1.8)  in  the  form 


(1.  20) 


For  a pure  homogeneous  plane  deformation,  referred  to  its  principal 


-12- 


(1.21) 


axes,  (1.  1)  become 

ya^^a^a  sum), 

in  which  the  coefficients  X.  are  positive  constants  and  are  readily  identified 

Of 

as  the  associated  principal  stretches.  If  the  material  is  isotropic,  the  co- 
ordinate axes  — in  this  particular  instance  — are  at  the  same  time  principal 
axes  common  to  the  now  constant  symmetric  tensor  fields  C,  G,  and  as 
well  as  to  a,  which  at  present  is  also  symmetric.  On  adopting  the  notation 


w w 


(1.  22) 


whenever  W is  to  be  regarded  as  a function  of  and  ^2  . one  finds  that 
(1 . 21)  reduces  (1.17),  (1.18)  to 


sum),  (a/^P), 


T =T  =T — (no  sum),  T o-O  (a^p), 

O'  oa  " »P 


(1.  23) 
(1. 24) 


where  T are  the  principal  Cauchy  stresses.  As  is  easily  seen,  (1.  24)  con- 
tinue  to  hold  locally  for  an  arbitrary  plane  deformation  (1.  1),  provided 
and  T are  the  local  principal  stretches  and  principal  Cauchy  stresses.  In 

Oi 

what  follows  we  will  need  to  refer  to  the  two-dimensional  Baker-Ericksen 
inequality,  according  to  which 

(Ti-T2)(Xi-X2)>0  if  ^1^>^2 

for  all  pure  homogeneous  deformations.  In  view  of  (1.  24),  this  condition  is 
equivalent  to 

(XiWi-X2W2)(Xi  -X2)>0  if  Xj^X^.  (1.  26) 
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We  shall  from  here  on  assume  that  the  strain-energy  density  and  all 
stresses  vanish  in  the  undeformed  state  and  thus  require  that  W(Xp\2)  obey 

W(l,  1)=0  . W^(l,  1)^0.  (1.27) 

The  transition  to  the  infinitesimal  theory  of  elastostatic  plane  strain,  which 
aims  at  a neighborhood  of  the  undeformed  state,  is  effected  through  a sys- 
tematic linearization  with  respect  to  the  displacement  gradients  of  the  finite 
theory  recalled  above.  The  underlying  limit  process,  which  corresponds  to 
X^"*l , confirms  that  and  coincide  to  dominant  order  in  this  approxi- 
mation. Upon  setting 

Cl=j[Wii(l,l)-Wi2(l,l)]  , X=Wj2(l,l),  (1.28) 

one  recovers  in  this  manner  from  (1.  17)  or  (1.  18)  and  (1.  7),  (1.  20)  the 
stress-displacement  relations  and  the  displacement  equations  of  equilibrium 
of  the  linear  theory  for  isotropic  materials  with  jj  as  shear  modulus  and  \ as 
Lame's  modulus.  Accordingly,  Poisson's  ratio  for  infinitesimal  deforma- 
tions is  given  by 

r Wj2(i,i) 

^"2Tx+kT)"w^j(l.l)+Wj2<^'^)’ 

With  reference  to  the  definitions  (1.  10)  and  (1.  12)  we  now  cite  from 
[2]  necessary  and  sufficient  conditions,  in  terms  of  restrictions  upon 
W(Xj,X2)  tbe  principal  stretches  X^,  for  the  ordinary  and  strong  ellip- 
ticity  of  the  two-dimensional  displacement  equations  of  equilibrium  (1.  7)  in 
the  presence  of  material  isotropy.  For  this  purpose  let 
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A=Xj-\2  . D = -(XjWj-\2W2)  if 


^ ‘^12^  “^2“^’ 


(1.  30) 


and  define  a symmetric  matrix  E=E(X^,X2)  as  follows:  for  all  X^>  0 such 

that  Xj/"X2. 


E =DW  (no  sum) 
aa  aa  ' 


1 W -W2  Wj+W2 

^12-^21  = 2 W^lWzz-  ^12  + T;TX^  ^12-T7Tx7 


(1  31) 


for  all  X >0  such  that  X,  =X-,=X, 
O l c 


E „=DW,  , . 

ap  11 


(1.32) 


Bearing  in  mind  the  symmetry  of  W(X^,  X2)  and  taking  for  granted  that  this 
function  is  twice  continuously  differentiable,  one  sees  that  D(Xj,  X2)  and 
£^^1»^2^  defined  by  (1.  31),  (1.  32)  are  continuous  on  the  entire  first  quad- 
rant of  the  principal- stretch  plane. 

As  sho\Mi  in  Section  2 of  [2],  for  an  isotropic  material  the  system 
(1.  7)  is  elliptic  at  a solution  and  at  a particular  material  point  if  and  only  if 
the  corresponding  local  principal  stretches  satisfy  the  inequalities 


Ej^jE22^0  , tiE  j + VE  I I E ^ > 0 , 


12  '^1 1 22 


(1.  33) 


where 


T|  = sgnE^^  (no  sum); 


(1.  34) 


further,  necessary  and  sufficient  for  the  strong  ellipticity  of  this  system  is 
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D>0  , Eji>0  . E22>0  , 


(1.  35) 


Note  from  (1,  30)  that  D>0,  if  is  equivalent  to  the  Baker-Ericksen  in- 

equality (1.26),  the  latter  being  a necessary  condition  for  strong  ellipticity. 
Moreover,  (1.  35)  evidently  imply  (1.  33),  as  should  be  the  case. 

At  the  undeformed  state,  characterized  by  strong -ellip- 

ticity conditions  (1.  35)  by  virtue  of  (1.  32),  (1.  30),  and  (1.  27)  reduce  to 


W^^(1,1)>0  . Wjj(l,  1)-Wj2(M)>0. 


(1.  36) 


These  inequalities,  in  turn,  because  of  (1.  28),  become 


|li>  0 , ^+2|!i>  0, 


(1.  37) 


which  are  precisely  the  familiar  conditions  for  the  strong  ellipticity  of  the 
linearized  displacement  equations  of  equilibrium  in  case  of  isotropy.  We 
proved  at  the  end  of  Section  2 in  [2]  that  when  (1.  37)  hold  true  — so  that 
strong  ellipticity  prevails  at  infinitesimal  deformations  — every  open  con- 
nected set  of  ordinary  ellipticity  in  the  (X^  that  contains  the  un- 

deformed state  (1,1)  is  necessarily  also  a domain  of  strong  ellipticity. 

Finally,  suppose  ellipticity  has  failed,  so  that  (1.  33)  are  violated  at 
some  material  point,  for  a particular  plane  deformation.  Then  there  must 
exist  at  least  one  (real)  characteristic  curve  passing  through  this  point. 

We  call  such  a curve  a "material  characteristic"  when  it  is  referred  to  the 
undeformed  configuration  and  apply  the  term  "spatial  characteristic"  to  the 
deformation  image  of  a material  characteristic.  Let  y be  the  angle  of  in- 
clination of  a local  spatial  characteristic  relative  to  the  first  principal  axis 
of  the  local  Cauchy  stress  tensor  T and  hence  also  of  the  left  deformation 


I 
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tensor  G.  As  shown  in  Section  3 of  [2],  cosZy  is  a solution  of 

^^1^1 1'’’^2^22‘^^1^2^12  2y 

+2(\2E22->^iEii)cos2y+(\^Ej^+\2E22+2>^1>^2^12)  <*• 

where  E(X.j,\2)  symmetric  matrix  introduced  in  (1.  31),  (1.  32).  The 

discriminant  of  this  quadratic  equation  is  non-negative  throughout  the  com- 
plement V of  the  domain  of  ellipticity  G,  characterized  by  (1.  33),  with  re- 
spect to  the  open  first  quadrant  of  the  (X^  ^ interior  of  V 

equation  (1.  38)  in  general  has  two  distinct  real  roots  within  the  interval 
[-1,1],  corresponding  to  four  distinct  spatial  characteristic  directions. 

The  latter  evidently  occur  in  two  pairs,  each  symmetrically  situated  with 
respect  to  y=0.  In  contrast,  for  points  on  the  common  boundary  of  C and  U 
(ellipticity  boundary),  which  are  associated  with  an  "incipient  failure  of  el- 
lipticity", there  is  at  most  one  such  pair  of  characteristic  directions,  the 

2 

two  real  roots  of  (1.  38)  being  necessarily  coalescent. 

2,  Two  special  elastic  materials  and  some  of  their  response  properties. 

For  future  illustrative  purposes  we  turn  here  to  a class  of  homoge- 
neous and  isotropic  elastic  materials,  whose  in-plane  response  to  a plane 
deformation  is  governed  by  an  elastic  potential  of  the  form 

W(I,J)  = ^[If(J)+g(J)]  (4>0),  (2.1) 

^ The  set  V of  non-elliptic  points  in  the  principal- stretch  plane  may  of 
course  be  empty  for  a particular  elastic  material. 

2 

See  the  exhaustive  discussion  of  incipient  failures  of  ellipticity  at  the  end 
of  Section  3 in  [2].  In  this  degenerate  case  there  may  exist  only  a single 
characteristic. 


LSfcr 


in  which  ^ is  a material  constant.  Note  that  this  potential  depends  linearly 
upon  the  invariant  I.  Since  the  corresponding  W(Xj,X.2)  is  to  satisfy  (1.27), 
the  functions  f and  g are  subject  to  the  requirements 

2f(l)+g(l)=0  . 2f(l)+2f'(l)+g'(l)=0.  (2.2) 

the  primes  denoting  differentiation.  For  every  fixed  choice  of  the  functions 
f and  g consistent  with  (2.  2)  one  can  readily  exhibit  an  infinite  class  of 
three-dimensional  elastic  potentials,  each  of  which  yields  (2.  1)  upon  spe- 
cialization to  plane  strain. 

We  now  particularize  (2.  1)  in  two  different  ways  and  hereafter  refer 
to  the  corresponding  materials  as  Material  1 and  Material  2.  Let 

f(J)=J  ^ , g(J)=2J-4  for  Material  1, 
f(J)=J  , g(J)=8j'^-10  for  Material  2, 

Accordingly,  for  Material  1, 

W(Xi,  X2)=^[2Xi\2+>^^^+X2^-4], 
whereas  for  Material  2, 

Both  materials  conform  to  (2.  2),  whence  the  energy  density  and  all  stresses 
vanish  in  the  undeformed  state.  Further,  the  first  of  (1.28)  gives  |a=|J  in 
either  instance,  so  that  u is  the  shear  modulus  for  infinitesimal  deforma- 
tions; on  the  other  hand,  the  second  of  (1.  28),  together  with  (1.  29),  furnish 

• 1 *2 

X =2^,  v = ^ for  Material  1 and  X =4p,  v = -g- for  Material  2.  Therefore,  (1.37) 
assure  strong  ellipticity  at  infinitesimal  deformations  for  both  materials. 


I (2.  3) 


(2.4) 


(2.  5) 


In  contrast,  as  will  become  apparent  in  what  follows,  neither  material  re- 
mains elliptic  at  all  plane  deformations. 

Our  next  objective  is  a detailed  discussion  of  relevant  response 
properties  pertaining  to  Material  1.  In  this  connection  we  note  first  that 
W(Xi,X.2)  (2.4)  is  the  restriction  to  plane  strain  of  an  elastic  potential 

proposed  by  Blatz  and  Ko  [3]  in  an  attempt  to  match  experimental  data  ob- 
tained in  tests  of  a highly  compressible  rubberlike  material.  The  idealized 
material  thus  adopted  in  [3]  was  the  subject  of  [1],  which  contains  a compre- 
hensive three-dimensional  treatment  of  its  response  characteristics  and  el- 
lipticity  restrictions.  Since  [1]  also  deals  explicitly  with  the  special  case 


of  plane  deformations  to  which  we  confine  our  attention  at  present,  the  sub- 
sequent results  are  mostly  cited  from  [1]  without  intermediate  detail. 

For  a pure  homogeneous  pi ane  deformation  of  the  form  (1 . 21)  one 
infers  from  (2.4)  and  (1.23),  (1.24)  the  in-plane  response 


* i 

appropriate  to  Material  1.  Specializing  (2.  6)  for  the  case  of  isotropic  plane 
. strain,  one  has 

- _4  -3 

Xj=X2=>^  . Tj=T2=t=u(1-X  ) , =02=a=u(X-X  ).  (2.7) 

•» 

t In  the  particular  instance  of  plane- strain  uni-axial  stress,  parallel  to  the 

Xj-axis,  (2,6)  furnish 

’ '^i=M(l-X[®)  , Qj  =u(Xj^-Xj^)  . (2.8) 

’ i 

We  consider  now  a homogeneous  plane  deformation  corresponding  to  a state 

I of  simple  shear  of  the  form 

3?,! 
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Yj^^Xj+Hx^  , Y2~^2’  (2.9) 

in  which  k is  a constant,  tan’^H  being  the  angle  of  shear.  The  principal 

stretches  of  this  deformation  are  related  to  h through  j 

! 

^1  ^ 

The  response  of  Material  1 to  such  a simple  shear  is  readily  deduced  with 
the  aid  of  (1.  17)  and  (1.  18).  In  this  manner  one  arrives  at 

^11"®’  ’ '’^12“''^21 

aii=a22=“U>t  , » a2^=u{x.+K  ). 


4.  L , \ 1 


As  is  clear  from  (2.  7),  in  isotropic  plane  strain,  t(X)  and  ct(X)  are 
concave,  monotone  increasing  functions  for  0<\<oo  and  T(X)-» -cd,  a(X)-*-oo  as 
\-»0,  but  T(X)-*u,  a(X)“*QO  as  X'^oo.  The  response  of  Material  1 to  plane-  strain 
uni-axial  stress  is  depicted  in  Figure  1.  According  to  (2.  8),  here  (X^ ) is  a 
concave,  monotone  increasing  function  for  0<X^<oo  and  Tj(Xj)-»-oo  as  X^-»0, 
while  T^(Xj)-»uas  X^'^oo.  In  contrast,  the  axial  nominal  stress  aj(Xj)  in- 
creases  steadily  from  aj^{0+)=-oo  to  a positive  maximum  at  Xj^=3'*-2.28  and 
thereafter  diminishes  steadily  toward  zero.  The  behavior  of  Material  1 in 
simple  shear  is  immediate  from  (2.  11).  Both  the  actual  shear  stress 
and  the  nominal  shear  stress  a^2  directly  proportional  to  the  amount  of 
shear  h,  but  is  a convex  strictly  increasing  function  for  0<k<cd.  A 

significant  nonlinear  effect  is  reflected  in  the  fact  that  the  shear  deforma- 
tion (2.  9)  induces  an  actual  normal  stress  T22I  the  latter  is  compressive 
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for  the  material  under  consideration. 

F rom  (2.  4),  (1.  30),  (1.31),  as  well  as  the  continuity  of  D(Xj , \^) 
^ap^^l'^2^'  fo*"  Material  1 at  all  , 


D = 


2 

^1^2 


' ^11=70 


,2  - 2 
= 3i_>0  , 


A 


M^2 


22"  2,  6 
^1^2 


>0 


M2-  12) 


^12”^,  6, 6 ^®^1^2‘^1'^2^ 


2X^X2 


J 


Hence  T|=1  in  (1.  34)  and  the  first  of  the  conditions  of  ordinary  ellipticity 
(1.  33)  holds  true  for  all  plane  deformations  of  this  material,  as  does  the 
Baker-Ericksen  inequality  (1.  26).  On  the  other  hand,  (2.  12),  the  second 
of  (1.  33),  and  (1.  14)  reveal  that  ellipticity  prevails  at  present  if  auid  only  if 


2 2 4 4 T 

14XiX2“^1 -X2>0  or  y<4. 


(2.  13) 


either  of  which  is  equivalent  to 


(2.  14) 


Moreover,  in  view  of  (1.  35),  these  conditions  are  also  necessary  and  suffi- 
cient for  strong  ellipticity. 

Let  (r,  9)  be  polar  coordinates  in  the  (X^,  X2)-plane,  defined  by 


X2=rcos0  , X2=rsin0  (0<r<oo,  O<0< ir/2). 


(2.  15) 


Then  (2.  14)  may  be  written  as 


0<r<oo  , . 


(2.  16) 


This  wedge-shaped  domain  of  ellipticity  C in  the  principal-stretch  plane  is 


, . \y  .♦ 


< 
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shown  in  Figure  2.  The  figure  also  displays  the  "isotropic  extension  path" 

1 

^1~^2'  "uni-axial  stress  path"  and  the  "simple- shear  path" 

supplied  by  (2.  7),  (2.  8),  and  (2.  10).  Since  the  first  of  these  paths 
lies  wholly  in  6,  Material  1 does  not  suffer  a loss  of  ellipticity  in  isotropic 
extension.  On  the  other  hand,  each  of  the  other  two  paths  intersects  the 
boundary  of  G twice.  Ellipticity  is  seen  to  fail  in  both  uni -axial  tension  and 
compression  according  as 

X^sp-*i2.  68  or  37;  (2.17) 

the  limits  of  ellipticity  for  uni-axial  stress  are  also  marked  in  Figure  1.  Fi- 
nally, in  simple  shear  a breakdown  in  ellipticity  is  found  to  occur  whenever 
either  principal  stretch  equals  or  exceeds  the  value  p ^=1.  93,  which  corre- 
sponds to  a shear  angle  of  approximately  55*. 

According  to  (1.  38)  and  (2.  12),  the  inclination  of  the  spatial  charac- 
teristics relative  to  the  first  principal  axis  of  Cauchy  stress,  for  Material  1 
obeys 

CO 

When  ellipticity  has  failed,  so  that  (2.  13)  is  violated,  (2.  18)  evidently  fur- 
nishes two  distinct  pairs  of  real  characteristic  directions,  except  on  the 
boundary  of  G,  where  Y = ±ir/4  are  the  only  solutions  of  (2.  18)  within 
Thus,  at  an  incipient  failure  of  ellipticity  the  spatial  characteristic  direc- 
tions for  Material  1 coincide  with  the  local  directions  of  the  lines  of  maxi- 
mum and  minimum  actual  shearing  stress.  Because  of  (2. 15),  equation  (2.18) 
alternatively  becomes 


X1-X2 


r 


cos2y  = 


/l-4sin  29 

cos20 


(2.  19) 


from  which  one  infers  that  the  characteristic  directions  do  not  vary  along 
any  ray  0=constant  in  the  non-elliptic  part  of  the  (\ j , \2)-plane.  This  com- 
pletes the  discussion  of  Material  1. 

A parallel  discussion  of  Material  2,  the  plane-strain  elastic  poten- 
tial of  which  is  given  by  (2.  5),  is  analytically  more  awkward.  Here  (2.  6) 


give  way  to 


Hence  for  isotropic  plane  strain. 


Xl=X2=X  , T^=T2=T=2^(X^-X■^)  . Oj  =a2=a=2u(X^-X’^).  (2. 

In  the  case  of  plane -strain  uni -axial  stress,  parallel  to  the  x^-axis,  one 
has  by  virtue  of  the  first  of  (2,  20), 


T2=a2=0  , 3X2+X^-4(XjX2,)'®  =0, 


the  last  of  which  implicitly  determines  the  transverse  stretch  as  a func- 
tion of  the  eixial  stretch  X^-  Thus, 

X2=X(Xj)  (0<Xj<oo)  (2.23 

and  (2.  20),  (2.  22),  (2.  23)  lead  to 


Tj=u[Xj-X^(Xj)]  , Oj=n[XjI(Xj)-X^(Xj)] 
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A 


Although  is  not  obtainable  in  elementary  form,  one  shows  without  dif- 

ficulty that  this  function  is  strictly  decreasing.  Further,  it  is  easy  to  de- 
duce an  explicit  parametric  representation,  in  terms  of  the  polar  angle  0 
introduced  in  (2.  15),  for  the  curve  corresponding  to  (2,  23).  This  repre- 


sentation was  used  to  plot  the  uni-axial  tension  path  appearing  in  Figure  4, 
as  well  as  the  stress- stretch  curves  depicted  in  Figure  3.  For  a simple 
shear  deformation,  characterized  by  (2.9),  (2.  10),  one  finds  at  present. 


'^11-^^  ' "^22  = ^ • '^12''^21=^^’^* 
2 3 


The  preceding  equations  reveal  certain  essential  qualitative  differ- 
ences between  Material  2 and  Material  1,  as  far  as  their  response  to  the 
special  homogeneous  deformations  under  consideration  is  concerned.  The 
functions  t(\)  and  a(\)  in  (2.  21),  which  govern  the  response  of  Material  2 
to  isotropic  plane  strain,  are  both  monotone  increasing  and,  as  is  the  case 
for  Material  1,  T(X)“» -oo,  o(X)"* -oo  for  X~»0;  however,  now  both  the  actual  and 
the  nominal  stress  tend  to  infinity  as  X"^co.  Similarly,  it  is  clear  from 
Figure  3 , which  pertains  to  uni-axial  stress,  that  the  true  axial  stress 
Ti(Xi)"»oo  as  Xj"»oo,  while  the  behavior  of  the  nominal  stress  Oj(X^)  is  quali- 
tatively the  same  as  in  Figure  1.  Further,  we  observe  on  the  basis  of  (2.25) 
that  a simple  shear  of  Material  2 induces  non-zero  actual  stresses  and 
^22’  which  are  tensile. 

In  place  of  (2.  12)  one  obtains  for  Material  2 at  all  X^^O, 
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D=hX^\2>0, 


2 2 2 *2 
E =3^  X,X,+(X,\.)^X‘ 

aa  1 2 ' 1 2'  a 


(no  sum). 


(2.  26) 


In  view  of  (2.  26),  (1.  30),  the  Baker-Ericksen  inequality  (1. 26)  is  satisfied 
also  for  Material  2.  Moreover,  an  appeal  to  (1.  33),  (1.34  ),  (1.  35)  con- 
firms that  ordinary  and  strong  eilipticity  once  again  prevail  if  and  only  if 
the  last  of  (1.  35)  holds  true.  The  latter  may,  with  the  aid  of  (2.  26)  and 
(2.  15),  be  written  as 


(10u)^-1)C+3+6ujVuu^C^+C+1>0, 


(2.27) 


provided  one  sets 


1 5 ~ 

UJ  = ysin20  , Q=r  uj® . 


(2.  28) 


Upon  exclusion  of  the  extraneous  root  of  the  quadratic  equation  in  Q obtained 
by  squaring  (2.  27),  one  arrives  at  a necessary  and  sufficient  condition  of  ei- 


lipticity for  Material  2 in  the  form 


(1  -2sin20)(sin20)* 


(o<0<^  , j2'<e<  j). 


(2.29) 


where  (r,  0)  are  the  polar  coordinates  in  the  principal- stretch  plane,  de- 
fined by  (2.  1 5). 

Figure  4,  which  is  the  counterpart  for  Material  2 of  Figure  2,  shows 
the  domain  of  eilipticity  £ appropriate  to  the  second  material.  This  figure 
also  exhibits  the  deformation  paths  in  the  (XjA2)"P^®^^®  corresponding  to 
isotropic  extension,  simple  shear,  and  uni-axial  stress  — the  first  two  of 
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which  are,  of  course,  the  same  as  in  Figure  2.  Here  again  ellipticity  never 
fails  in  isotropic  extension  and  breaks  down  for  all  sufficiently  severe  sim- 
ple siiears.  In  contrast  to  the  results  for  Material  1,  however,  those  given 
in  Figure  4 reveal  the  analytically  ve rifiable  fact  that  Material  2 cannot  suffer 
a loss  of  ellipticity  in  uni-axial  compres sion;  the  failure  of  ellipticity  in  uni- 
axial tension  is  seen  to  occur  at  an  axial  stretch  \^-3,  04  (see  also  Figure  3). 

The  inclination  y of  the  real  spatial  characteristic s for  Material  2, 
whose  existence  is  assured  if  the  ellipticity  condition  (2.  29)  is  violated, 
follows  from  (1.  38)  by- recourse  to  (2.  26)  and  (2.  15).  In  this  manner  one  is 
led  to 


cos2y 


1 

4(|cos29 


1^  2-2 
-16uu  C > 


(2.  30) 


0)  and  Q being  the  auxiliary  functions  of  the  polar  coordinates  (r,  0)  adopted 
in  (2.  28).  Equation  (2.  30)  is  the  analogue  for  Material  2 of  (2.  19).  At  all 
points  in  the  complement  of  £ with  respect  to  the  open  first  quadrant  of  the 
(Xl , X^)  "P^^3.ne , except  for  the  points  on  the  boundary  of  £,  (2.  30)  yield  four 
distinct  characteristic  directions.  On  the  boundary  of  £ there  results  a sin- 
gle pair  of  distinct  characteristics,  whose  inclincations  are  determined  by 


cos2y 


2 - sin 2 0 
2cos20  ’ 


(2.  31) 


where  0 is  the  angular  polar  coordinate  of  the  boundary  point  in  question. 
Thus,  the  characteristic  directions  at  an  incipient  failure  of  ellipticity  of 
Material  2,  unlike  those  associated  with  Material  1,  are  no  longer  constant. 
In  particular,  at  an  incipient  breakdown  of  ellipticity  in  uni-axial  tension 
parallel  to  the  Xj-axis,  one  has  0^3.  51  (see  Figure  4)  and  (2.31)  yields 
Y ^ ± 9.  46  . 


J 


T 
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3.  Elastostatic  fields  with  discontinuous  deformation  gradients.  Piece - 
wise  homogeneous  elastostatic  shocks. 

At  this  stage  we  turn  to  our  main  objective,  which  concerns  the  ex- 
istence and  nature  of  a class  of  elastostatic  fields  with  discontinuous  defor- 
mation gradients  appropriate  to  a homogeneous,  but  not  necessarily  iso- 
tropic, hyperelastic  solid.  Suppose  such  a solid,  in  its  undeformed  config- 
uration, occupies  the  entire  (Xj , X2)-plane  R spanned  by  a rectangular  car- 
tesian coordinate  frame  X.  Let  X,  with  the  unit  direction  vector  j^,  be  a 
straight  line  through  the  origin  of  X,  so  that 

£:2S=x(§)=L?  (-oo<?<oo),  (3.1) 


and  call  N the  unit  normal  vector  of  X obtained  by  a counter-clockwise  ro- 

+ 

tation  of  L through  a right  angle.  Next,  we  designate  by  R and  R the  two 

open  half-planes  into  which  X divides  R,  with  the  understanding  that  points 
+ 

into  R(  Figure  5)  Consider  now  a piecewise  homogeneous  plane  deformation 
of  the  form 


+ + + + 
for  all  X in  R , J =detF>0 

I 

Fx  for  all  x in  R , j =detF>0  . 


(3.  2) 


+ 

Here  F and  F are  constant  (nonsingular)  tensors,  which  evidently  represent 

+ 

the  position  - independent  deformation-gradient  fields  prevailing  on  R and  R, 
respectively,  while  J and  J are  the  corresponding  Jacobian  determinants. 

According  to  ( 1 . 6),  the  nominal  stress  field  induced  by  the  deforma- 
tion (3.  2)  is  given  by 


+ 3-  + _ - _ 

a=Wp(F)  on  R , on  R 


ti 


mis 


(3.  3) 
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and  automatically  satisfies  the  equilibrium  equation  (1.  5)  on  either  side  of 
X.  We  shall  assume  that  the  displacement  field  associated  with  (3.  2)  is 
continuous  across  X,  so  that 

for  every  x on  X,  (3.  4) 

which  — in  view  of  (3.  1)  — is  equivalent  to 

FL=FJL.  (3.  5) 

This  assumption  permits  us  to  extend  the  mapping  (3.2)  continuously  onto 
X and  excludes  from  our  present  considerations  any  separation  or  gliding 
of  material  along  the  unique  deformation  image  X^^  of  X,  Moreover,  as  is 
clear  from  (3,  1),  (3.  3), 

(-oo<§<oo).  (3.6) 

Figure  5(b)  illustrates  typical  deformation  images  of  the  three  rectangles 
shown  in  Figure  5(a)  under  the  mapping  (3.2),  subject  to  the  continuity  re- 
quirement (3.  4). 

Finally,  equilibrium,  i.  e.  the  balance  of  forces  across  X,  demands 
the  continuity  of  the  Piola  tractions  at  X,  so  that 

aN=oN5s,  (3.7) 

where  s is  the  nominal  traction  exerted  on  the  material  in  R by  the  mate- 
+ 

rial  occupying  R in  the  undeformed  configuration. 

When  F^F,  we  shall  refer  to  the  elastostatic  field  characterized  by 
(3.  2),  (3.  3),  together  with  the  continuity  conditions  (3.  5),  (3.  7),  as  a piece- 
wise  homogeneous  elastostatic  shock  (equilibrium  shock);  further,  we  shall 
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call  the  straight  lines  £ and  governed  by  (3.  1)  and  (3.  6),  the  material 
and  the  spatial  shock-line. 

The  first  question  that  arises  in  connection  with  piecewise  homoge- 
neous elastostatic  shocks  as  characterized  above  concerns  their  existence: 

are  there  such  shocks  for  a given  hyperelastic  material  and  a given  defor- 

+ 

mation  gradient  F?  If  so,  how  many?  Is  their  existence  contingent  upon  } 

restrictions  on  the  governing  elastic  potential?  Further,  in  case  such  I 

shocks  exist,  what  are  the  corresponding  orientations  of  the  shock-lines  j 

1 

and  values  of  F?  Also,  what  is  the  nature  of  the  emerging  elastostatic  i 

fields  and  what  kind  of  field  discontinuities  at  the  shock-lines  do  they  entail?  5 

t 

Before  we  can  attempt  to  find  at  least  partial  answers  to  the  fore- 
going questions  we  need  to  explore  in  some  detail  the  kinematics  of  equi- 
librium shocks.  To  this  end  we  designate  by  i the  unit  direction  vector  of 
note  on  the  basis  of  (3.  6)  that 

f = cFL=cFL  , c=-i  - =— , (3.8)  i 

IfTT  lEfel 

I 

and  assign  to  the  unit  normal  vector  n resulting  from  a counter-clock- 
wise rotation  of  I through  ti/2.  Clearly,  separates  the  two  open  half- 
planes ftj^and  that  are  the  deformation  images  of  ^ and  R;  also,  n points 
+ 

into  R^(see  Figure  5).  Now  let  $ and  cp  stand  for  the  angles  of  inclination, 
relative  to  the  x,  -axis,  of  £ and  £^,  respectively.  Both  of  these  angles, 
hereafter  referred  to  as  the  material  and  spatial  shock-angles,  may  be 
confined  to  the  interval  [ -it/ 2 , u / 2]  and 

Lj  =cos4>  , L2=sin4>,  Nj  = -sin$,  N^^cos^*, 
jf^=coscp  , f2=sincp  , n|^  = -sincp  , n2=coscp. 


(3.  9) 
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L =&•  oNq  , N Lo  , = e =60  (3.  10) 

a ap  p ' Of  Pof  p Of  QfP  p O'  Por  P 


Further,  frorri  (3.  8),  (3.  10)  cind  the  last  of  (I.  16)  follows 


1+-1  I'-l  1+T  1-T 

L = -(F)  "f=-(F)  , N=^F  n-^F  n. 


(3.  11) 


Because  of  (1.  13),  the  left  deformation  tensors  of  the  homogeneous  defor- 

+ 

mations  on  R and  ft  are  given  by 


+ ++T  - --T 

5=EE  , G=FF^, 


{3.  12) 


whence,  squaring  the  last  of  (3.  11),  one  has 


2 n • Gn  n • Gn 


(3.  13) 


Also,  (3.  8)  and  (3.  9)  yield  the  relation 


^21 

Fii+Fj^ta„* 


(3.  14) 


between  the  material  and  the  spatial  shock-angle. 

Since  £ and  F are,  by  hypothesis,  nonsingular  tensors  with  positive 
determinants,  there  is  a tensor  B such  that 


-4-1  - + 

B=F(F)  , F=BF  , detB>0. 


(3.  15) 


One  may  therefore  resolve  the  mapping  (3.  2)  into  the  two  successive  plane 
deformations: 


z=Fx  for  all  x in  ft. 


(3.  16) 
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z for  all  z in 


Bz  for  all  z in  R, 


(3.  17) 


We  shall  call  the  mapping  (3.  16),  which  is  homogeneous  on  the  entire  plane 
Si,  the  intermediate  deformation  belonging  to  the  original  deformation  (3.  2) 
and  refer  to  (3.  17)  as  the  supplementary  deformation.  The  displacement 
continuity  condition  (3.  4)  evidently  demands  that  the  supplementary  defor- 


mation carry  into  itself,  i.  e. 


(3.  18) 


as  is  also  apparent  from  (3.  5)  and  (3.  15),  (3.  11).  At  this  stage  we  intro- 
duce a second  coordinate  frame  with  the  same  origin  and  the  unit  base 
vectors  (f  ,n),  so  that  X'  is  obtained  by  a rotation  of  the  original  frame  X 
through  the  spatial  shock-angle  cp  (see  Figure  5).  On  referring  (3.  18)  to 
the  frame  X'  one  finds  at  once  that  the  displacement  field  of  the  deforma- 
tion (3.  2)  is  continuous  if  and  only  if 


B'^(  B'  ) = 


(3.  19) 


where  K and  6 are  two  as  yet  arbitrary  constants,  whose  kinematic  signifi- 
cance will  emerge  presently.  Meanwhile  we  take  note  of  the  matrix  relation 


B=R^E!'r  , G'=RGR  , R = 


cos!p  smcp 


-sincp  coscp 


(3.  20) 


^Boldface  letters  carrying  a prime  stand  exclusively  for  the  matrix  of 
scalar  components  in  X' of  the  corresponding  tensor  or  vector. 
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and  thus  conclude  with  the  aid  of  (3.  9)  that  (3.  19)  is  equivalent  to 


B fo+Ki  no  + 6n  n„. 
QfP  o p cr  p a p 


(3.21) 


Furthermore,  from  (3.  Z),  (3.15),  (3.  ZO)  follows 


+ + 

J =detF  =detBdetF=JdetB^ , 


(3.ZZ) 


so  that  by  virtue  of  (3.  19), 


- + + 
J‘*6J  f 


(3.  Z3) 


1 


if  principal  stretches  of  the  deformation  on  R and  hence  also  of 

the  intermediate  deformation  (3.  16).  Next,  (3.  IZ),  (3.  15)  give 


G=BGg'^  , G'=B'g'(B')^, 


(3.  24) 


the  second  of  which,  in  conjunction  with  (3.  19),  leads  to 


^1 1 ~^1 1 ‘'’^^^IZ''’^  ^ZZ  ’ ^ZZ”^  ^ZZ’ 


^1Z“^Z1  ^^IZ’*’^^  ZZ^  • 


y (3.25) 


As  for  the  physical  interpretation  of  the  parameters  h and  6,  we  in- 
fer first  from  (3,  15),  (3.  19),  (3.  ZO)  that 


F=^  , B=B'=^  if  and  only  if  H=0  , 6=1, 


(3.  26) 


in  which  case  the  supplementary  deformation  (3.  17)  is  the  identity  mapping 
and  the  original  deformation  (3.  Z)  is  trivial  in  the  sense  of  being 


r 


Recall  (1.  14). 


i • 
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homogeneous  on  ft,  and  thus  no  longer  constitutes  an  equilibrium  shock. 


Moreover,  (3.23)  and  the  principle  of  mass  conservation  imply 

44.., 


(3.  27) 


where  p and  p are  the  mass  densities  of  the  material  occupying  ft^j,  and  ft^j^  in 
the  deformed  configuration.  Thus  (3.  27)  supplies  an  interpretation  of  6 in 
terms  of  the  area-ratios  or  the  ratio  of  mass  densities  appropriate  to  an 
equilibrium  shock. 

An  additional  kinematic  meaning  of  6,  as  well  as  a geometric  inter- 
pretation of  the  parameter  k,  come  into  evidence  if  one  factors  the  matrix 
of  (3.  19)  as  follows: 


1 H 1 r 1 0|  r 1 H 


0 6j  LO  6 10  1 


(3.  28) 


This  factorization  reveals  that  the  supplementary  deformation  (3.  17)  ad- 
mits a decomposition  into  a simple  shear  of  amount  h,  parallel  to  the  z ^ - 
axis  (spatial  shock-line),  followed  by  a uni-axial  stretch,  with  the  stretch- 
ratio  6,  at  right  angles  to  this  axis.  The  preceding  resolution  is  illustrated 
in  .Figure  6:  the  diagram  in  Figure  6(b)  depicts  the  supplementary-deforma- 
tion history  of  the  two  unit  squares  shown  in  Figure  6(a),  which  pertains  to 
the  intermediate  configuration.  It  should  be  emphasized  that  Figure  6 is 
based  on  6<1,h>0  and  requires  obvious  modifications  if  6^1  or  k^O. 

For  future  purposes  we  mention  here  the  important  special  case  in 

+ 

which  the  deformation  on  ft  in  (3.  2)  is  pure  homogeneous,  having  X as  a 
principal  frame.  IfX.j,X^  are  once  again  the  corresponding  principal 
stretches,  one  has  in  this  particular  instance 
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+ 

F: 


’^1 

0 

+ 

0 

7 

0 

1 

0 

^2J 

(3.29) 


while  (3.  14)  reduces  to 


tancp  = -r — tan4>. 
^1 


(3.  30) 


We  note  also  that  the  second  of  (3.  29)  holds  true  for  an  arbitrary  homoge- 
neous deformation  on  provided  X is  a principal  frame  for  its  deforma- 
+ 

tion  tensor  G. 

We  return  now  to  the  traction  continuity  condition  (3.  7),  which  — on 
account  of  (3.  3)  — is  equivalent  to 


W„  (F)N=W„  (F)N. 


(3.  31) 


Suppose  T and  _T  represent  the  constant  Cauchy  (actual)  stress  field  to  which 

+ 

the  piecewise  homogeneous  deformation  (3.  2)  gives  rise  on  and  re- 
spectively. Then,  by  (1. 4)  and  (3.  3), 


+ 1 ,+  . + T + .1  ,-,-T 

(F)F  on  , T = — Wj,(f)F  on 
T ~ J ~ 


(3.  32) 


An  appeal  to  (3.  1 1 ),  the  validity  of  which  depends  on  the  continuity  of  the  dis- 
placements, at  once  confirms  that  (3,  31),  (3.  32)  imply 


(3.  33) 


l! 


where  ^ denotes  the  Cauchy  (actual)  traction  exerted  by  the  material  occu- 
pying R,j,  in  the  deformed  configuration  on  the  material  in  Conversely, 

(3.33),  (3.  32),  (3.  1 1)  assure  that  (3.31),  and  hence  (1.  7) , holds  true. 
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Thus,  in  the  presence  of  continuous  displacements,  the  nominal  tractions 
are  continuous  ac  ross  the  material  shock -line  £ and  only  ^ the  actual 
tractions  are  continuous  across  the  spatial  shock -line 

For  later  convenience  we  cite  at  this  place  also  the  particular  form 
assumed  by  the  traction  continuity  condition  (3.  33)  in  the  event  the  hyper- 
elastic material  under  consideration  happens  to  be  isotropic.  In  this  case 
(3.  33),  because  of  (1.  18),  furnishes 


|Wj(I,  J)  G+Wj(I,  J)  n=  4Wj(i  , j)  5+Wj(i,j)_i  n,  (3.34) 

T J 


+ + - _ + - 
provided  I,  J and  I,  J stand  for  the  scalar  invariants  of  G and  G,  so  that  in 

accordance  with  (1.  14), 


+ +22^  I ^ 

I=trG=Xj+\2  > J=vdetG  =detF=X^\2> 


I=trG  , J =v  detG  =detF. 


(3.  35) 


It  is  clear  from  the  second  of  (3.  15),  together  with  (3.  9)  and  (3.  21), 


F-F(cp,  H,  6)=B(cp,  K,  6)F. 


(3.  36) 


Hence,  bearing  in  mind  the  developments  leading  up  to  (3.  21),  one  sees 

that  the  displacement  continuity  condition  (3.  5)  alone  constrains  the  possi- 

+ 

ble  values  of  F,  in  a shock  corresponding  to  a fixed  prescribed  F,  to  a 
three-parameter  family.  The  existence  of  such  a shock  therefore  hinges 
on  the  existence  of  a spatial  shock-angle  cp  in  [ -it/2,  tt/2],  a real  value  of 
the  shear  parameter  H,  and  a positive  value  of  the  stretch  parameter  6, 
such  that  the  traction  continuity  condition  (3.  33)  has  a solution  F(cp,K,  6) 
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other  than  the  trivial  solution  supplied  by 

- - + + 

F=F(cp,0,  l)=B(cp,  0,1)F=F  (--rr/Z^cpSir/Z).  (3.37) 


Furthermore,  since  (3.  33)  constitute  only  two  scalar  restrictions  on  the 
three  parameters  (cp,  H,  6),  one  would  anticipate  that  U there  exists  an  equi- 
librium shock  for  a given  deformation  gradient  F,  there  exists  a one-param- 
eter family  of  such  shocks. 

We  now  prove  the  following  theorem,  which  establishes  a necessary 
condition  for  the  existence  of  shocks  of  the  kind  under  consideration: 

If  there  exists  ^ piecewise  homogeneous  elastostatic  shock  in  a hyper- 
elastic material,  then  the  displacement  equations  of  equilibrium  associated 
with  this  material  must  suffer  a los s of  strong  elliptic ity  at  some  homogene- 
ous deformation. 

+ 

With  a view  toward  establishing  this  claim,  let  F and  F be  the  defor- 
mation gradients  of  the  existing  shock,  so  that  from  (3.  2),  (3.  5),  and  (3.  31), 
+ _ + _ + 

detF>0  , detF>0  , FL=FL  , W^(  F)  N=W„(  F)  N,  (3.38) 


where  L and  N are  the  unit  direction  and  the  unit  normal  vector  of  the  mate- 
rial shock-line  £.  Next  define  a family  of  tensors  by  means  of 


F(a)=aF+(l-a)F  (O^cySl). 


(3.  39) 


Then, 

E(0)=F  , F(1)=F,  (3.40) 


while  (3.  39)  and  (3.  15)  give 

■4 

F(a)^[(H+(l -a)B]F  (O^asI).  (3.41) 


L 
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From  this  relation,  in  turn,  because  of  (3.20),  (3.  19),  (3.27),  and  the 
first  of  (3.  38),  follows 

+ 

detF(o')  =(or6  + l -a)detF>0  (Osa^l),  (3.42) 

whence  F(a)  is  admissible  as  a deformation-gradient  tensor  of  a homogene- 
ous plane  deformation  for  every  a in  the  interval  [0,  1]. 

According  to  the  third  of  (3.  38),  there  is  a vector  K such  that 

Let  0 be  the  scalar -valued  function  defined  by 

+ 

0(a)=K.[Wj,(F(a))-Wp(F)]  N (O^a^l),  (3.44) 


so  that  owing  to  (3.  38)  and  (3.  40), 


0(0)  =0(1)  =0. 


(3.45) 


Thus,  by  virtue  of  the  mean-value  theorem,  there  is  a number  a in  (0,  1) 
such  that 


d0 

da 


a=a 


=0. 


(3.46) 


But  (3.  44),  (3,  46)  in  conjunction  with  (3,  39),  (3.  43),  and  (1.  8)  lead  to 
Finally,  set 

F=F(S)  , M--K/ IkI  (3.48) 


and  appeal  to  (1.9)  to  arrive  at 
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M.Q(N;F)M=0.  (3.49) 

which  contradicts  the  strong -ellipticity  condition  (1.  ll).  This  completes 
the  proof. 

It  remains  a matter  of  conjecture  as  to  whether  or  not  the  existence 
of  an  equilibrium  shock  necessitates  not  only  a loss  of  strong  ellipticity  but 
also  a failure  of  ordinary  ellipticity,  as  is  intuitively  plausible.  In  this 
connection  we  refer  once  more  to  a result  obtained  in  [2]  and  cited  in  Sec- 
tion 2 of  the  present  paper.  ^ It  follows  from  this  result  that  for  an  isotropic 
hyperelastic  material  a failure  of  strong  ellipticity  implies  a failure  of  or- 
dinary ellipticity,  if  the  set  of  all  points  in  the  principal-stretch  plane  at 
which  strong  ellipticity  prevails  is  a domain,  i.  e.  open  and  connected,  and 
includes  the  undeformed  state.  ^ 

A related  comment  pertains  to  the  existence  of  a piecewise  homoge- 
neous equilibrium  shock  within  the  linear  theory  of  homogeneous  and  iso- 
tropic elastic  solids.  One  verifies  easily  that  in  this  setting  such  a shock 
exists  if  and  only  if  the  elastic  constants  satisfy 

u/^0  and  \+2|a=0  or  v = l,  (3.  50) 

which  require  a failure  of  ordinary  ellipticity  of  the  linearized  displace- 
ment equations  of  equilibrium. 

4.  Weak  piecewise  homogeneous  elastostatic  shocks. 

Taking  for  granted  the  existence  of  piecewise  homogeneous  equilib- 
rium shocks  in  the  (possibly  anisotropic)  hyperelastic  material  under 

^See  the  discussion  following  (1.  37), 

We  recall  from  Section  2 that  these  conditions  are  met  for  the  two  special 
isotropic  materials  discussed  there.  See  also  Figure  2 and  Figure  4. 


• 
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[ consideration,  we  confine  our  attention  in  this  section  to  shocks  that  are 

+ 

I weak  in  the  sense  that  F remains  close  to  F.  Indeed,  motivated  by  (3.  36), 

i (3.  37),  and  the  observations  following  (3.  37),  we  assume  here  that  there 

! 

exists  a one-parameter  family  of  shocks  corresponding  to  the  given  gradi- 

+ 

ent  F,  depending  on  the  parameter 

e=l-6,  (4.1) 

and  suitably  smooth  near  e=0. 

Specifically,  we  presuppose  that  there  are  functions  cp(e)  and  H(e), 
both  twice  continuously  differentiable  in  a neighborhood  of  e=0,  such  that 
F(e)  defined  by 

+ + 

F(e)=B(q)(e).H(e),  6(e))F=B(e)F,  (4.  2) 

where 

6(e)  = l-e  , K(0)=0  (4.3) 

and  B(cp,  H,  6)  is  the  tensor  characterized  by  (3.  21),  conforms  to  the  traction 

i 

I continuity  condition  (3.  33)  — or,  equivalently,  to  (3.  31)  — throughout  the 

I 

I 

1 neighborhood  at  hand.  Since  at  present  the  trivial  (shockless)  solution  (3.37) 

^ corresponds  to  e =0  and  e is  evidently  a measure  of  the  departure  from  this 

► 

solution,  we  shall  henceforth  refer  to  e as  the  shock- strength  parameter. 

The  kinematic  significance  of  e in  terms  of  the  supplementary  deformation 
is  immediate  from  that  of  6 (see  Figure  6);  also,  (4.  1)  and  (3.  27)  yield 
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existence  of  the  e -family  of  shocks  postulated  above,  to  ihe  lowest  signifi- 
cant order  in  €.  For  this  purpose  we  observe  first  that  the  material  shock- 
line and  the  spatial  shock-line  are  now  e -dependent.  We  thus  write  Z{e) , 
in  place  of  and  consequently  also  L(e) , N(e ) , i (e ) ,£(e ) , as  well  as  ft(e), 

R(e) , R,,.(e) , Byje)  and  4>(e).  Furthermore,  any  function  previously  defined  on  R 
or  ft.,,  and  dependent  on  the  value  of  the  gradient  F is  at  present  to  be  regarded 
as  a function  of  e.  With  this  understanding  all  of  the  results  in  Section  3 up 
through  (3.  35)  hold  true  identically  in  e. 

In  view  of  the  assumed  smoothness  of  cp(e),K(e)  and  by  virtue  of  the 
second  of  (4.  3),  one  has  the  Taylor  expansions 


N(e)=N(0)+N(0)e+o(e)  , n(e)^n(0)+n(0)e+o(e), 
K(e)  =h(0)g+o(6)  as  e-»0; 


y (4.  5) 


here  and  in  what  follows  a dot  placed  above  a letter  indicates  differentiation 
with  respect  to  the  parameter  e.  On  the  other  hand,  using  (3.  Zl),  (3.9),  as 
well  as  (4.  3),  one  finds  that  B(e)  obeys 


provided  £ is  the  vector  defined  by 

q-H(0)^(0)-n(0).  (4.7) 

F rom  (4,  6)  follows 

and,  keeping  (4.  Z)  in  mind,  one  thus  draws 


(4.  9) 
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Consequently,  an  appropriate  two-term  Taylor  expansion  leads  to 

aw(F(e))  aw(F)  a^W(F) 

^^ny  (4.10) 

Finally,  entering  (3.  31)  with  (4.  10)  and  using  the  first  of  (4.  5),  as  well  as 

(3.  11),  one  sees  that  the  traction  continuity  condition  is  fulfilled  up  to  the 
order  of  e if  and  only  if 

2 + 
d W(F) 

(4.11) 

This  equation,  in  turn,  may  be  written  as 

Q(N(0);F)^.0, 

where  Q(N;F)  is  the  symmetric  tensor  introduced  through  (1.8),  (1.  9). 
Now^/tO  according  to  (4.  7).  Therefore  (4.  12)  implies  that 

detQ(N(0);F)=0. 

Drawing  on  the  discussion  of  (1.  10),  (1.  11)  in  Section  1,  one  is  thus  led  to 
the  following  conclusions; 

(i)  A necessary  condition  for  t^  existence  of  a one -parameter  fam- 
iliC  of  e^uili^^  (of  the  kind  under  present  consideration)  .s  that 

jjsplacemen_t  equations  of  equilibrium  associated  with  the  hyperelastic 

— suf^  a loss  of  ordinary  ellipticity  ^ ^ given  homogeneous  de- 
formation on  ft. 

(n)  In  ^ wea^- shock  i.  e,  ^ e^O,  the  material  shock-line 

J(e)  ai^^  spatial  shock -W  j!^(c),  respectively,  tend  to  a material  and 
a spatial  characteristic  associated  with  the  homogeneous  defo rmati on  on  ft. 


1 
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In  compa ring  conclusion  (i)  with  the  theorem  proved  at  the  end  of 

Section  3,  we  note  that  the  present  result  yields  a breakdown  of  ordinary 

ellipticity  — rather  than  merely  of  strong  ellipticity;  on  the  other  hand  (i) 

presupposes  the  existence  of  an  entire  family  of  elastostatic  shocks,  while 

in  the  theorem  referred  to  but  a single  such  shock  was  required  to  exist. 

Conclusion  (ii)  restricts  the  number  of  possible  distinct  limiting 

material  (and  spatial)  shock-lines  to  the  number  of  distinct  characteristic 

directions  admitted  by  (1.  11),  which  cannot  exceed  four.  In  particular,  if 

the  material  is  isotropic  and  X is  chosen  to  be  a principal  frame  for  £ and 

T,  the  inclination  cp(0)  of  the  limiting  spatial  shock-line  Z^(0)  must  be  such 

that  cos2cp(0)  is  a real  root  in  the  interval  [-1,1]  of  the  quadratic  equation 

(1.  38)  for  cos2y  , with  y replaced  by  cp(0).  As  an  illustrative  example  con- 

+ 

sider  the  isotropic  Material  1 discussed  in  Section  2 and  suppose  that  F in- 
duces an  incipient  failure  of  ellipticity.  Then,  as  is  clear  from  the  remarks 

following  (2.  18),  cp(0)=±tT/4  so  that  £,,.(0)  necessarily  coincides  with  a trajec- 

+ 

tory  of  the  maximum  actual  shear  stress  of  the  constant  stress  field  ^ pre- 

+ 

vailing  on  R. 

A kinematic  interpretation  of  the  vector  q introduced  in  (4.  7)  is 
readily  arrived  at.  Let  w(^;g)  be  the  family  of  supplementary  displacements 
corresponding  to  the  e -family  of  equilibrium  shocks  at  hand.  Then,  by  (3.17) 
and  (4.  8), 

+ 

w(^;g)=0_  on  ft^:.(e),  (4.14) 

w(^;e ) = [B(g) = fn(0) +o(g)  on  lRjj,(G).  (4.15) 

From  (4.  15)  and  tlic  second  of  (4.  5)  one  infers  that 
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(4.  16) 


provided  ^w(^;e)  is  the  tensor  field  with  the  components  9w^(^;e)/9Zp  in  the 

frame  X.  For  future  convenience  we  now  adopt  the  following  notation:  if 

+ - + - 

f is  a scalar  or  a vector  field  defined  on  R and  R (or  on  R^  and  R^)  that  suf- 
fers a finite  jump  discontinuity  across  £ (or  across  £*),  we  write 


[f];=f-f 


(4.  17) 


+ - + - 

for  the  jump  in  f as  £ (or  £j,J  is  traversed  from  R to  R (or  from  R^j^  to  R^|,). 
Equations  (4,  16),  (4.  14)  thus  justify  the  assertions 


[V;^]^=^e+o(e)  , [7w^]'=0. 


(4.  18) 


Hence  q^e  is  the  lowest-order  approximation  to  the  jump  across  £^(e)  in 
the  derivative  of  w^(^;e)  normal  to  Z^{e),  whereas  the  directional  deriva- 
tive of  w^(z^;e)  parallel  to  Is  continuous  at 

+ 

While  the  characteristic  lines  of  the  deformation  on  R that  emerge 
upon  a loss  of  ordinary  ellipticity  are  associated  with  jumps  in  the  second 
normal  derivatives  of  the  displacements^  equilibrium  shocks  are  seen,  on 
the  basis  of  (4.  18),  to  involve  discontinuities  in  the  first  normal  derivatives 
of  the  displacements  across  the  shock-line.  Also,  it  should  be  emphasized 
that  £,|,(e)  is  in  general  not  a characteristic  line  of  the  deformation  on  R(e), 
although  it  tends  to  such  a line  as  e-»0,  i.  e.  in  the  weak -shock  limit. 

According  to  (4.  1),  (4.  2),  and  (3.  23),  the  jumps  across  £(e)  in  the 
area-ratio  and  the  mass  density  obey 

[J]'  = -je  = -\j\2e  , [p]'=pe.  (4.19) 

^Recall  the  role  of  the  characteristics  reviewed  in  connection  with  (1.  11). 
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In  order  to  compute  the  jump  across  £(e)  in  the  strain-energy  density,  we 
first  rely  on  (4.  9)  to  see  that 


+ aW(J') 

W(F(e))  =W{  F)  + 3P  ~ q^n^(0)F^pg  + o(e). 


(4.20) 


On  account  of  (1.  4)  and  (1.  6),  this  equation  may  be  written  as 

+ + + 

W(F(e))  =W(F)  +Jjg.j;n(0)e  + o(e) 


(4.21) 


or,  by  virtue  of  (3.  33)  and  (4.  17), 


_ + + 

[W]“  =J;^(0)  *qe+o(e)  , t(0)=Tn(0). 


(4.22) 


We  calculate  next  the  jump  across  S.^(e)  in  the  Cauchy  stresses  T^. 
From  (3.3  2),  (4.  19)  follows 

- , 1 aw(F(e))-  - + 

' J(e)=(l-e)J.  (4.23) 


These  formulas,  together  with  (4.  9),  (4.  10),  and  (1.  8),  permit  one  to  de- 
duce 


Finally,  it  is  of  interest  to  compute  the  jump  across  £^(6)  of  the 
scalar  normal  stress  acting  parallel  to  the  spatial  shock-line,  i.  e.  of 
if  T^p  are  the  components  of  ^when  the  latter  is  decomposed  in  the  frame  X' 
(see  Figure  5).  Clearly, 


(4.  25) 


r 
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and  (3.  10)  yield 


i^(e)ip(e)=i^(0)ip(0)  + o(l)=6^p-n^(0)np(0)  + o(l). 


(4. 26) 


Substitution  from  (4.  24)  and  (4.  26)  into  (4.  25),  after  a lengthy  calculation 
that  makes  use  of  (4.  7),  (3.  11),  (1.  9)  and  the  traction  continuity  condition 
in  the  form  (4.  12),  eventually  leads  to  the  result: 


[T'li]  '=q  •^(0)e+o(e), 


(4.27) 


where  uu=(ju(F),  while  cu(F)  is  the  tensor  with  the  components 


(4.28) 


in  the  frame  X and  JT(F)  the  actual  stress  tensor  associated  with  F through 
(1.4),  (1.6). 

The  lowest-order  jump  approximations  (4.  18),  (4.22),  (4.24),  and 

+ 

(4.  27)  involve  — beyond  quantities  fully  determinable  from  the  given  F and 

the  known  response  function  W(F)  — the  vector  originally  introduced  in 

(4.  7).  Furthermore,  this  vector  involves,  in  addition  to  the  unit  direction 

vector  f (0)  and  the  unit  normal  vector  n(0)  of  the  spatial  characteristic 

£(0),  also  the  still  unknown  value  h(0). 

We  now  determine  A.(0)  from  (4.  12),  and  for  this  purpose  exclude 

+ 

the  degenerate  case  in  which  Q (N(0);F)  is  the  null  tensor.  Thus  assuming 


QhQ(N(0);F)/^0, 


(4.29) 


we  show  first  that 


k(0)=0  if  and  only  if  Qn(0)=^. 


(4.  30) 
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n(0)*Q^(0);t0  if  Qn(0);t0, 


(4.  31) 


Now  (4,  7),  (4.  12),  because  of  (4.  30),  (4.  31),  furnish 


. . ja(O)  *Qs,(0) 

H(0)=0  ifQn(0)-0,  if  Qn(0)9^0. 


(4.  32) 


where  Q again  abbreviates  Moreover,  (4.  32)  and  the  last  of  (4.5) 


n(0)*Qn(0) 


(4.  33) 


as  a lowest-order  weak-shock  approximation  to  the  amount  of  shear  inherent 
in  the  supplementary  deformation. 

Our  next  task  is  to  specialize  some  of  the  foregoing  results  for  the 
case  of  material  isotropy.  In  this  instance  one  finds  with  the  aid  of  (1.9) 
and  (1.  20)  that,  for  every  unit  vector  N, 


Q (N;F)=2W,6^  + 
ay ' I ay 


4W„F  «N«F  „N„+J^W^.rF«^NoF‘^N  + 
II  ap  (3  YP  p JJ  p«  p pY  p 


2jw,^(f:^n„f  n +f:^n„f  N ). 
IJ'  pof  p YP  P Py  P “P  P 


(4.  34) 


To  see  this,  suppose  h(0)=0.  Then  Qjn(O)  =£  follows  at  once  from  (4.7), 

(4.  12).  Next,  if  k(0)f^0,  Qn(0)  cannot  vanish  either  since  then  ^ and  n(0) 
would  be  linearly  independent  null  vectors  of  the  tensor  Q,  so  that  Q would 
have  to  be  the  null  tensor,  contrary  to  (4.  29).  Hence  (4.  30)  is  true.  One 
confirms  similarly  that 


We  now  use  (4.  34)  in  (3.  66),  invoke  (3.11),  (3.12),  (3.13),  as  well  as  (3.9), 
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+ 

and  choose  X as  a principal  frame  for  G,  so  that  the  component  matrix  of 
+ 

G in  X is  given  by  the  second  of  (3.  29).  In  this  manner  we  eventually  ar- 

+ 

rive  at  the  following  result  for  h(0)>  provided  (4.  29)  holds  and  Q(N(0);F)^(0) 
also  fails  to  vanish: 

2+2  + 

2aW  +4a  W +J^W  .+4aJW 

k(0)= i (4.35) 

2(2aW^+JW^j)b 

where 

a=X^sin^cp(0)  + X.2Cos^cp(0)  , b-(\  j -\2)sincp(0)coscp(0)  (4.  36) 

and  the  partial  derivatives  of  W with  respect  to  the  deformation  invariants 

+ + 

are  understood  to  be  evaluated  at  (I,  J),  the  latter  being  supplied  by  (3.  35). 

If  one  adheres  to  the  above  special  choice  of  the  frame  X,  (4.  22)  is 
in  the  present  circumstances  found  to  imply 

[W]‘ = - (2dWj+JWj)e+o(e),  (4.  37) 

provided 

+2  + 

2aW  +J  W +2aJW 

d-a  + OH{0)--- ^ (4.38) 

2(2aWjj+JWjj) 

and  a is  given  by  the  first  of  (4.  .ih).  The  tensor  introduced  in  (4.  28)  in 
the  case  of  the  isotropic  hyperelastic  solid  turns  out  to  have  the  components 

= (4.  39) 

Finally,  (4.  27)  and  (4.  39)  enable  one  to  confirm  that 


•■'Or  ' 
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4 


I T'l  J - T 

j 


jj+lJW^j]e+o(€),  (4.40) 


if  the  previous  choice  of  X is  retained.  Also,  the  derivatives  of  W appear- 

+ + 

ing  in  (4.  37),  (4.  38),  (4.  40)  stand  for  their  corresponding  values  at  (I,  J). 

For  future  purposes  we  record  here  the  weak-shock  approximations 
of  K(e)  for  the  two  special  isotropic  materials  discussed  in  Section  2.  Re- 
calling (2.  1),  (2.  3),  we  note  that  the  strain-energy  densities  governing  the 
reponse  of  these  materials  obey: 


'y 

W(I,  J)  = ^(IJ'  +2J-4)  (u>0)  for  Material  1, 

W(I,  J)  = ^{IJ+8J‘^-10)  (u>0)  for  Material  2. 

On  the  basis  of  (4.  35),  (4.  41)  and  (4.  5)  one  finds  that 

K(e)  = 2b  (I-a-)e+o(e)  for  Material  1, 

K(e)=  - — (a+J’®  )e+o(e)  for  Material  2, 


y (4.41) 


J 


1 (4.  42) 


with  a and  b given  by  (4.  36). 

In  connection  with  the  weak- shock  jump  estimates  deduced  in  this 

section,  it  is  essential  to  recognize  that  the  signs  of  some  of  these  jumps 

may  depend  not  only  on  the  particular  material  considered,  but  also  on  the 

+ 

particular  nature  of  the  pre-assigned  deformation  on  ft;  moreover,  the  sign 
of  each  jump  depends  on  the  sign  of  the  shock- strength  parameter  e.  The 
sign  of  e,  in  turn,  cannot  be  determined  in  the  absence  of  information  be- 
yond that  contained  in  our  present  characterization  of  piecewise  homoge- 
neous elastostatic  shocks.  In  Section  6 we  shall,  on  energetic  grounds. 
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propose  an  additional  requirement  that  leads  to  a removal  of  this 
indeterminacy. 

5.  Equilibrium  shocks  of  finite  strength  for  a particular  isotropic  hyper- 
elastic material. 

We  turn  now  to  an  instructive  illustrative  example  concerning  the 
global  existence  and  character  of  piecewise  homogeneous  elastostatic 
shocks  in  a special  (homogeneous)  isotropic  hyperelastic  material.  The 
following  analysis  is  based  on  Material  1 of  Section  2,  the  strain  - ene  rgy 
density  of  which  is  given  by  the  first  of  (4.  41).  To  simplify  this  analysis 
we  shall  assume  here  from  the  start  that  the  second  of  (3.  29)  is  in  force, 

so  that  X is  a principal  frame  for  the  deformation  tensor  G of  the  given  de- 

+ + 
formation  on  and  hence  for  the  actual  stress  tensor  ^ as  well. 

The  traction  continuity  condition  (3.  34)  at  present  reduces  to 

lt(G-Il)n  = -i(G-il)n.  (5.1) 

*3  - j3 

Upon  referring  (5.  1)  to  the  frame  X'  introduced  in  Section  3 (see  Figure  5) 
and  bearing  the  first  of  (3.  35)  in  mind,  one  arrives  at 


On  setting 

(q>0.X2>0),  (5.3) 

one  obtains  by  means  of  the  second  of  (3.  29)  and  the  second  of  (3.  20): 
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G'ii=i(Xi+\2)^(l-2pcos2cp+(3^), 


G'22=^(Xl+X2)^(l+2pcos2cp+p^), 


+ / 1 2 

Gi2  = ^(Xi+X2)  Psin2cp. 


(5.  4) 


The  components  G^^  are  now  computable  from  (3.  25),  while  J is  related  to 
J through  the  first  of  (3.  27).  In  this  manner  (5.  2)  lead  to  the  two  scalar 
traction  continuity  conditions 


2 6^Psin2cp  1 +2pcos2q)+p^)  +2Psin2cp, 


6^(1  -2Pcos2cp+p^)  = 


1 -2Pcos2cp+p^+K^(  1 +2pcos2cp+P^)+4Hpsin2cp. 


(5.  5) 


Observe  that  the  given  deformation  on  51,  in  the  example  under  con- 
sideration, enters  (5.  5)  exclusively  through  the  principal  - stretch  param- 
eter p adopted  in  (5.  3).  Also,  (5.  3)  and  (3.  27)  require  that  P and  6 satisfy 


•1<P<;1  , 6=1 -e>0. 


wr.ere  e is  the  shock  - strength  parameter  introduced  in  (4.  1).  From  here 
on  we  shall  give  preference  to  e over  6.  The  trivial  solution  of  (5.  5),  which 
signifies  the  absence  of  a shock,  is  furnished  by  e =0  (6=1),  H=0,  and 
-ir/2scp^iT/2.  We  now  seek  the  answer  to  the  following  question.  For  what 
values  of  P in  (-1,1)  and  e<l  do  there  exist  nontrivial  solutions  H=K(e;P)  and 
cp=cp(e;P)  of  the  simultaneous  equations  (5.  5),  such  that  H is  real  and  cp  in  the 
interval  [ -ir/2  , ir  / 2]?  Moreover,  we  wish  to  find  all  such  solutions. 

It  is  clear  from  the  structure  of  (5.  5)  that  if  a pair  of  values  (H,cp) 


J ‘ 
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satLsfLes  these  two  equations,  then  so  does  the  pair  Consequently, 

the  limitation  of  the  spatial  shock-angle  to  the  range 

O^cp^TT/Z  (5.7) 

entails  no  loss  in  generality.  From  the  first  of  (5.  5)  and  the  second  of 
(5.  6)  one  has 

^^2(6.^-l)PsinZy 
1 +2(3cos2(p+p 

and  substitution  for  k from  (5.  8)  into  the  second  of  (5.  5),  in  view  of  (5.  6) 
and  (5.  7),  is  found  to  yield 


sin2cp  = - 


1-1 

P 

2|p 

|k(e) 

k(e): 


r 2 

6+6  + 1 


>0 


6=1 -e. 


(5.  9) 


if  once  and  for  all  we  remove  the  indeterminacy  in  cp  at  e =0  by  requiring 
continuity.  Since  q)  is  to  be  real,  the  positive  right  member^  of  the  first  of 
(5.  9)  cannot  exceed  unity.  Hence, 


^k(e)  (-oo<e<l)  , 1p1<1. 


(5.  10) 


These  two  inequalities  are  easily  seen  to  be  equivalent  to 

h(e)^lP|<l  . h(e)  = [k(e)+yk^(e)  + l (-oo<e<l)  (5.11) 


and  are  necessary  in  order  that  (5.  5)  admit  a solution  (K,cp),  subject  to 
(5. 6),  (5.  7),  for  a given  pair  of  values  (e=l-6,P).  Conversely,  if  (5.  11)  is 
fulfilled,  such  a solution  is  supplied  by  (5.  8),  (5.  9). 


1 


Recall  (5.  6). 


-51- 


Thus,  in  the  present  circumstances,  an  elastostatic  shock  of 

strength  e<l  exists  if  and  only  if  the  principal  stretches  A.j,X.2  inherent  in 

+ 

the  pre-assigned  deformation  on  ft,  which  determine  the  parameter  p 
through  (5.  3),  are  such  that  the  point  with  the  rectangular  cartesian  coor- 
dinates (e,  1P|)  of  the  (e,  lp|)-plane  lies  in  the  "admissible  region"  G char- 
acterized by  (5.  11)  and  shown  in  Figure  7.  In  this  connection  we  note  on 
the  basis  of  the  definitions  of  k(e)  and  h(e)  in  (5.  9),  (5.  1 1)  that 


k 

o 


=k(0)=h(0)=-^, 

73 


A 


k(l-)=0  , h(l-)=l  , h(-oo)=0  , k(l-)=h(l-)=0. 


V (5.  12) 


k(e)<0  , h(e)>0  (-oo<e<l),  J 


provided  primes  denote  differentiation  with  respect  to  the  shock- strength 
parameter  e. 

According  to  (5.  9),  every  pair  of  values  (e,P)  admitted  by  (5.  10)  or 
(5.  11)  gives  rise  to  two  spatial  shock-angles  in  the  interval  [0,ir/2]: 

Tl(€;P)=|-sin 


2lPlk(e) 


(O^cPjSir/4) 


V (5.  13) 


The  amounts  of  (supplementary)  shear,  and  > associated  with  the  re- 
spective shock-angles  cpj  and  cp^  are  readily  found  by  recourse  to  (5.  8)  and 
(5.  13).  In  this  manner  one  obtains 

1 

H^(e:(3)  = p^(^)  (a  = l,2)  , 6 = 1-g,  (5.14) 

where 


liiiSiiMi 


/ * 
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R^(e;P)=sgnpi^k(e)+(-l)'^^^Q(e:P)  . Q(e;P)  = 
" l-p 


4p\^(e) 

— 2'T 

L (i-P  ) 


-1 


(5.  15) 


We  observe  that  Q(e:P)  is  real  whenever  (e,P)  conform  to  the  admissibility 
requirement  (5.  10)  and  that 

<Pj(e:P)  = cp^(e;p)  = j , Kj(€;P)=K2(e:p)  if  =k(e),  (5.16) 

i.  e.  , if  (e , I P 1 ) is  a point  on  the  boundary  of  the  admissibility  region  G dis- 
played in  Figure  7.  Further,  as  is  apparent  from  (5.  14),  (5.  15), 

«-j(e;-P)= -H2(€;P)  , >t2(e;-P)= -»<-i(e;P).  (5.17) 

With  reference  to  the  remark  preceding  (5.  7)  we  emphasize  that  (5.  13), 

(5.  14)  in  general  determine  four  distinct  equilibrium  shocks  appropriate  to 
the  particular  material  at  hand. 

One  confirms  easily  with  the  aid  of  (2.  13),  (3.  35),  and  (5.  4)  that  the 

displacement  equations  of  equilibrium  are  elliptic  at  the  homogeneous  de- 

+ 

formation  given  on  R if  and  only  if 

|p|<k  (5.  18) 

o ,/J 

+ 

For  the  sake  of  brevity  we  shall  say  that  "ellipticity  prevails  on  R”  when- 
ever (5.  18)  holds  true.  Evidently,  such  is  the  case  at  all  points  (e , 1p1)  of 
the  admissible  region  G (Figure  7)  that  lie  below  the  line  |pl=k^,  while  all 

t* 

points  of  G on  or  above  that  line  correspond  to  a failure  of  ellipticity  on  ft. 

In  particular,  such  a failure  of  ellipticity  on  ft  occurs  for  all  admissible 
pairs  (e  , I P I ) with  e^O. 

We  now  seek  to  ascertain  for  what  points  in  G the  ensuing  homoge- 
neous deformation  on  ft  is  elliplii-  or  otherwise. 


Aciording  to  (2.  13), 
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ellipticity  prevails  on  if  and  only  if 

i<4j. 


(5.  19) 


On  appealing  to  (3.  35),  (3.25),  (3.23),  (5.4),  remenibe ring  that  G and  G 
have  the  same  trace,  we  find  with  the  aid  of  (5.  8)  — after  a rather  lengthy 


computation  — that 


. (Xj+X^) 


2 r 


(l+(3^)(l+6^)+2P(6^-l)cos2cp 


4P^(6'^-l)sin^2cp 


l+2pcos2cp+(3 


, J -6X  j X^ > 


(5.20) 


in  which  cp  = cp  (e:P)  is  furnished  by  (5.  13).  Owing  to  (5.  20)  and  (5.  3),  the 


ellipticity  condition  (5.  19)  for  ft  becomes 

( 


1 


46(1-(3^)L 


1 +(3^ ) ( 1 + 6 ^ ) +2  p ( 6^  - 1 ) c o s 2cp+ 

l+2pcos2q)+p  J 


<0. 


(5.21) 


The  inequality  (5.21),  together  with  (5.6),  enables  one  to  prove  that  there 
i s ^ failure  of  ellipticity  on  ft  fo r all  admis sible  pairs  (e  , | P | ) with  e ^ 0 ; in 
contrast,  the  subregion  of  G corresponding  to  6>0  contains  points  (6,  |pl)  at 
which  ellipticity  prevails  on  ft,  as  well  as  points  at  which  it  does  not. 

For  reasons  that  will  become  apparent  later  on,  we  confine  our  at- 
tention in  the  remainder  of  this  section  to  shocks  of  strength  e^O. 

Figure  8,  which  relies  on  (5.  13),  (5.  14),  illustrates  the  dependence 
of  the  two  spatial  shock-angles  cp^(e;P)  within  the  interval  [0,tt/2],  and  of  the 


associated  amounts  of  shear  >t^(e:P),  upon  e for  e^O  and  P-0.  65.  At  this  val- 


ue of  P the  maximum  possible  e is  approximately  0.  22,  as  is  clear  from 


Figure  7.  Both  k,  and  h,  arc  negative  for  0<eSe  . Observe  that,  in 
° 12^  max 


The  schematic  inset  diagram  at  the  top  of  Figure  8 corresponds  to  a posi- 
tive  value  of  k. 


j 
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in  agreement  with  (5.  16),  cp,  =cp^  -it/4  and  K,  =k~,  when  e -e  , in  which  case 
(e,  1(3  1)  is  a point  on  the  boundary  of  the  admissibility  region  G.  While  cpj, 
cp^,  and  Hj  vary  monotonically  with  e,  the  curve  for  in  Figure  8 exhibits  a 
minimum  in  the  interior  of  the  interval  fO,e  1. 

The  jump  in  the  strain-energy  density  W,  across  the  appropriate  ma- 
terial shock-line,  is  deducible  from  the  first  of  (4,  41)  by  means  of  (3.  35) 
and  (5.  40).  The  result  of  this  computation  may  be  put  into  the  form 


L‘-p5i.H6).LU66 

6(l+2Pcos2cp+p  ) 


where  cp  = cp^(e;P)  is  available  from  (5.  13)  and  Of-1  or  a=2  depending  on  the 
particular  shock  at  hand.  Since  6=1 -e>0,  it  follows  that  all  shocks  of  posi- 
tive strength  e possible  in  Material  1 give  rise  to  a decrease  in  the  energy 

+ 

density  as  the  shock-line  £ is  traversed  from  R to  R. 

Because  of  (5.  12),  (5.  13), 


- 2 ' 

cp^(0;p)=-|sin~^ 

It  is  not  difficult  to  verify  by  recourse  to  (5.  3)  that  equations  (2.  18)  are 
satisfied  if  y equals  either  of  the  limiting  spatial  shock-angles  in  (5.  23). 
Hence,  in  the  limit  as  e"*0,  the  spatial  shock-lines  appropriate  to  Material  1 
tend  to  the  spatial  characteristics  associated  with  the  prescribed  homoge- 
neous deformation  on  R,  the  latter  being  accompanied  by  a loss  of  ellipticity. 
This  conclusion  reflects  a general  result  concerning  weak  elastostatic 
shocks^  established  in  Section  4.  Similarly,  a Taylor  expansion  about  e =0 
applied  to  (5.  14)  confirms  the  consistency  of  this  global  result  with  its  weak- 
shock  counterpart  in  (4.  42). 


See  the  discussion  of  (4.  13). 


-55- 


6 • Time-dependent  piecewise  homo)Jeneous  elastostatic  shocks.  Energy 
considerations . Dissipativity . 

The  loss  of  ellipticity  of  the  field  equations  of  elastostatics  at  par- 
ticular — sufficiently  severe  — deformations  of  certain  hyperelastic  solids  is 
analogous  to  the  change  of  type  that  may  occur  in  the  partial  differential  equa- 
tions governing  steady  irrotational  flows  of  a compressible,  inviscid  fluid^. 
These  equations  are  elliptic  or  hyperbolic  at  a point  of  the  flow  field  accord- 
ing as  the  corresponding  particle  velocity  is  subsonic  or  supersonic. 

One  of  the  important  features  of  compressible  flows  is  the  possible 
occurrence  of  shock- surfaces  across  which  the  fluid  pressure,  density  and 
velocity,  as  well  as  the  entropy,  suffer  jump  discontinuities.  The  simplest 
example  of  this  kind  in  the  theory  of  steady  plane  flows  is  that  in  which  a 
plane  of  discontinuity  separates  two  uniform  flows.  Further,  this  discontin- 
uous flow  field  is  a close  analogue  of  a piecewise  homogeneous  elastostatic 
shock.  The  mathematical  counterpart  in  the  foregoing  fluid-flow  problem  of 
the  conditions  of  displacement  and  traction  continuity  are  restrictions  aris- 
ing  from  the  balance  of  mass,  momentum,  and  energy  across  the  gas-dynam- 
ical  shock-plane.  The  ensuing  one-parameter  family  of  shocks  may  be  re- 
ferred to  the  shock- strength  parameter 


See  Courant  and  Friedrichs  [?]  for  a general  treatment  of  compressible 
flows.  A detailed  discussion  of  the  strict  cinalogy  between  steady  gas  flows 
and  anti-plane  shear  in  finite  elastostatics  may  be  found  in  [s]. 

2 

Actually  there  are  two  symmetrically  located  families  of  shocks;  in  contrast, 
the  higher-order  elasticity  problem  gives  rise  to  two  pairs  of  symmetrically 
situated  one-parameter  families  of  equilibriimi  shocks  — at  least  in  the  exam- 
ple of  Section  5. 
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where  p and  p are  the  respective  fluid  densities  on  the  upstream  and  down- 
stream side  of  the  shock.  Equation  (6.  1)  is  identical  with  the  second  of 

(4.4),  in  which  p and  p are  the  mass  densities  of  the  elastic  material  occupy- 
+ - 

ing  R..,  and  R,,,,  respectively,  in  the  deformed  state. 

In  gas  dynamics  the  shock  conditions  mentioned  above  are  accom- 
panied by  the  independent  requirement  that  the  entropy  of  a fluid  particle 
shall  increase  as  the  shock- surface  is  traversed.  This  entropy  condition 
models  the  dissipative  nature  of  the  process  of  shock-formation  in  the  ab- 
sence of  viscosity;  it  leads  directly  to  the  conclusion  that  g must  be  positive 
and  hence  p>  p . 

The  role  of  entropy-like  conditions  was  studied  extensively  by  Lax 
[9]^  for  quasi-linear  hyperbolic  systems  of  conservation  laws  in  two  inde- 
pendent variables.  Lax's  work  is  motivated  prirr  arily  by  the  initial  - value 
problem  for  a system  of  partial  differential  equations  in  which  time  appears 
explicitly  as  one  of  the  two  independent  variables. 

In  view  of  our  remarks  concerning  stationary  shocks  in  gas  dynam- 
ics, as  well  as  on  independent  physical  grounds,  it  is  natural  to  subject  elas- 
tostatic  shocks  likewise  to  an  additional  limitation  that  assures  their  dissi- 
pative character.  For  this  purpose  we  find  it  essential  to  generalize  the  no- 
tion of  an  elastostatic  shock  defined  in  Section  3 to  a time -dependent  family 
of  such  equilibrium  shocks.  It  should  be  made  clear  that  time  will  play  merely 
the  part  of  a history  parameter  in  this  context  since  there  are  no  inertia  ef- 
fects involved  in  the  following  quasi- static  considerations. 


^This  paper  contains  references  to  related  earlier  investigations  by  the  same 
author. 
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As  in  Section  3,  let  R stand  for  the  entire  (Xj , x^)-plane  spanned  by 
a fixed  (time  - independent)  rectangular  cartesian  coordinate  frame  X.  Let 
£{t)  (t^^t^t^)  be  a time  -dependent  family  of  straight  lines,  which  conforms 
to  the  parametrization 


<C(t) : x=x(5,  t)  =p^(t)  + 5L(t)  ( -oo<§<oo,  t j^t^t^). 


(6.2) 


Here  L(t)  is  the  orienting  unit  direction  vector  of  £(t)  at  the  time  t,  while 
p^(t)  is  the  instantaneous  position  vector  of  a point  that  is  attached  to  £ and 
participates  in  the  rigid  motion  of  £ relative  to  the  frame  X;  evidently,  % is 
the  directed  distance  from  this  point,  measured  along  £ (see  Figure  9).  We 
take  for  granted  that  £ and  L are  continuously  differentiable  functions  of  the 
time  on  the  interval  [tj,!^].  Further,  paralleling  the  agreements  intro- 
duced at  the  beginning  of  Section  3,  we  call  N(t)  the  unit  normal  vector  of 
£(t)  resulting  from  a counter-clockwise  rotation  of  L(t)  through  ir/2.  Also, 

we  denote  by  R(t)  and  R(t)  the  open  half-planes  into  which  R is  divided  by 

+ 

£(t)  at  the  instant  t,  with  the  proviso  that  N(t)  points  into  R(t). 

At  this  stage  we  define  a time-dependent,  piecewise  homogeneous, 
family  of  plane  deformations  through 

"+  + + 

F(t)x+b(t)  for  all  x in  k{t)  (t^^t^t^) 

F(t)x+b(t)  for  all  x in  K(t)  (t^^t^t^), 

+ + - - , , 

J =detF>0  , J =detF>0  on  [t^t^]. 


These  equations  are  th^  time -dependent  counterpart  of  (3.2).  Note  that 
the  translation  vectors  are  needed  here  since  £(t),  given  by  (6.  2),  may 
not  pass  through  the  origin  at  all  times.  Actually,  but  for  a lack  of  sym- 
metry, one  of  these  two  vectors  could  have  been  omitted. 
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I 


1 


+ “ 

The  (position-independent)  deformation  gradients  F,  F,  as  well  as  the 

+ _ 

translation  vectors  Jb,^,  are  required  to  possess  continuous  first  time-de- 
rivatives on  [t|,t2]. 

Instead  of  (3.  3)  we  now  have 

+ + + - 

a(t)-W^{F(t))  on  R(t)  , a(t)  =Wj,(F(t))  on  R(t)  (t^^tst^),  (6.4) 


where  a and  £ are  the  nominal  stress  fields  produced  by  the  time-depend- 

+ 

ent  deformation  (6.  3)  on  R and  R,  if  the  solid  occupying  R in  its  undeformed 
configuration  is  hyperelastic  and  W is  its  strain-energy  density.  The  dis- 
placement-continuity condition  (3.  4)  at  present  gives  way  to 
+ + _ 

F(t)^+b(t)  =F(t)x+b(t)  for  all  x on  £(t)  (tj^t^t^):  (6.  5) 


because  of  (6.  3),  this  requirement  is  met  if  and  only  if 

which  take  the  place  of  (3.  5).  Assuming  the  mapping  (6.  3)  to  have  been 
extended  continuously  onto  X(t),  we  let  £.,,(t)  stand  for  the  deformation  im- 
age of  £(t)  (tj^tst^). 

The  stress  field  (6.4)  clearly  satisfies  the  equilibrium  equation  (1.  5) 

+ 

on  R(t)  and  R(t)  at  each  instant;  but  the  balance  of  force  across  £(t)  demands 
the  continuity  of  the  nominal  tractions^  at  <£(t); 

(6.7) 

+ 

When  F(t)/:F(t)  for  tj^t^t^,  we  shall  refer  to  the  family  of  elasto static 


Owing  to  the  continuity  of  the  displacements,  the  balance  of  mass  on  either 
side  of  £(t)  assures  the  mass  balance  across  £(t). 


C » 
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fields  characterized  by  (6.  3),  (6.  4),  together  with  the  continuity  conditions 

(6.  6),  (6.  7),  as  a time  -dependent  (piecewise  homogeneous)  elastostatic 

shock  with  the  (moving)  material  and  spatial  shock-lines^  £(t)  and  £j|,(t) 

(t^^t^t^).  Observe  that  at  any  fixed  instant  t such  a shock  is  an  equilibrium 

shock  of  the  type  defined  in  Section  3,  if  ^(t)=^  and  the  translation  vectors 

+ 

^(t),Jb(t)  also  vanish. 

It  is  readily  seen  that  equations  (3.  8)  through  (3.  15)  hold  identically 
on  [tjjt^]  for  time -dependent  shocks  provided  i,n,4>,cp  and  all  field  quantities 
retain  their  previous  meaning  but  are  now  regarded  as  functions  of  time.  On 
the  other  hand,  at  any  fixed  instant  the  mapping  (6.  3)  admits  the  resolution 

+ + ^ 
^=F(t)x+b(t)  for  all  ^ in  ft,  j 


fz  for  all  z in 


y (6.  8) 


B(t)z-B(t)b(t)+b(t)  for  all  z in  R,|,(t), 


where  ft,,Jt)  and  ft.„(t)  are  the  instantaneous  deformation  images  of  ft(t)  and 
ft(t)  associated  with  (6.  3),  while  B obeys  the  first  of  (3.  15)  on  [tj,!^].  Argu- 
ments strictly  parallel  to  those  employed  in  Section  3 confirm  that  equations 
(3.  18)  through  (3.  37)  are  valid  on  [tj,!^]  in  the  present  circumstances. 

Note  that  the  frame  X.'  with  the  base  vectors  (f  ,n)  is  now  time -dependent, 
as  are  the  parameters  K,  6,  and  e,  whose  kinematic  significance  remains 
unaltered.  Figure  6 is  applicable  to  the  instantaneous  supplementary  de- 
formation inherent  in  the  decomposition  (6.  8),  except  that  the  spatial  shock - 
line  while  parallel  to  — need  not  be  coincident  with  the  y*j-axis. 

Consider  a time -dependent  elastostatic  shock  of  duration  [tj,t2]  and 


The  term  "material  shock-line"  is  somewhat  misleading  in  the  present 
setting  since  £(t),  being  occupied  by  different  particles  at  different  times, 
is  not  a material  line  . 


r 
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let  be  in  the  open  interval  (t^jt^).  With  a view  toward  arriving  at  the  no- 
tion of  dissipativity  of  the  given  shock  at  time  t^  we  adopt  the  following 
definition.  A domain  of  the  (x^x^) -plane  is  a test  domain  admissible  at 
time  t^,  if  & is  the  interior  of  a fixed,  smooth  and  convex,  simple  closed 
curve  C and  the  material  shock-line  £(4^)  intersects  iB.  Then,  evidently, 

JC(t)  intersects  j0  for  all  values  of  t in  a neighborhood  of  t^,  which  we  shall 
call  a "time-range  of  intersection  of  Further,  for  every  t in  such  a 
range  of  time,  £(t)  intersects  the  convex  boundary  C of  j6  in  two  distinct 
points  — say  with  the  position  vectors  x(§^(t),  t)  and  t),  where 

indicated  in  Figure  9.  Also,  £(t)  divides  the  stationary  do- 

4* 

main  j6  into  the  two  time -dependent  sub-domains  f)(t)  and  &{t),  which  are  the 

+ - 

intersections  of  & with  the  half-planes  ft(t)  and  R(t),  respectively  (see  Figure  9). 

At  this  stage  we  choose  a test  domain  & admissible  at  the  given  instant 
and  let  U(t)  stand  for  the  total  strain  energy  in  a slab  of  unit  thickness  — at 
any  time  t within  a time- range  of  intersection  of  jS  — stored  in  the  material 
that  occupies  & in  the  undeformed  configuration.  Accordingly,  we  have 
+ +-  - + + 

U(t)=W(t)A(t)  + W(t)A(t)  , W(t)=W(F(t))  , W(t)=W(F(t)),  (6.9) 

+ - + - 
where  A(t)  and  A(t)  are  evidently  the  respective  areas  of  fi(t)  and  fi(t).  Dif- 
ferentiation of  the  first  of  (6.  9),  in  view  of  (6.  4)  and  the  second  of  (6.  9), 
yields  the  time-rate  of  change  of  the  strain  energy  U in  the  form 

U*=WA* +WA* +a«  f:  A+5o  F:  A.  (6.10)^ 

Py  Py  Py  Py 

From  the  divergence  theorem  follows  at  once 

^Here  and  in  the  sequel  dots  used  as  superscripts  indicate  differentiation 
with  respect  to  the  time. 
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, 5„pA(t)4x^Vpd/.,  (6,  11) 

al(t)  ajfi(t) 

+ - + 

if  ,s  is  the  arc -length  along  the  boundaries  9i6(t),  9jB(t)  of  the  domains  jB(t), 

£(t)  and  v is  the  appropriate  unit  outward  normal  vector.  On  the  other  hand, 

it  is  readily  shown,  and  intuitively  obvious,  that 

^§2(t) 

A’(t)= -A‘(t)= -j  V(|,t).N(t)d§,  (6.  U) 

§l(t) 

provided 

V(5.t)-|^x(?.t),  (6.13) 

so  that^  V(?,t)  is  the  velocity,  relative  to  the  frame  X,  of  the  point  on  the 
moving  shock -line  Z(t)  with  the  position  vector  x(5,t). 

We  now  use  (6.  11),  (6.  12)  in  conjunction  with  (6.  10)  and  then  invoke 
(6.2),  the  continuity  conditions  (6.  5),  (6.  7),  as  well  as  the  overall  equilib- 
rium of  the  material  occupying  &(t)  and  fi(t).  In  this  manner  we  are  led  to 

U -(t)  =js(x,  t)  -vCx,  t)46+| i[W(t)]'V(?,  t)  -N(t)  + s(t)  .[v(x,  t)i;  }d?.  (6.14) 

9j6  ?l(t) 

Here  s(x,t)  is  the  Piola  traction  on  the  outer  side  of  9jS,v(x,  t)  is  the  parti- 


cle velocity  given  by 
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4 


+ 

appropriate  field  values  as  £(t)  is  traversed  from  ft(t)  to  ft(t);  thus,  in  par- 
ticular, 

[w(t)];=w(F(t))-w(F(t))  , [v(x,t)];=[F-(t)];i(§,t)+[b-(t)];.  (6.16) 

From  the  second  of  (6.  16),  together  with  (6.  2),  (6.  5),  and  (6.  13)  follows 


te.t)];=-[F(t)];V(?,t). 


(6. 17) 


Substituting  from  (6.  17)  into  (6.  14),  and  appealing  once  more  to  the  traction- 
continuity  condition  (6.  7),  one  finds  after  elementary  manipulations: 


provided 


?2(t) 

U’(t)  =j£(x,  t)  *v(x,  t)c^+j  [P(t)]^(t)  •V(?,  t)d§, 

?j(t) 

P(F)=W(F)UF'^a(^:)  , P(t)-P(F(t))  , P(t)=P(F(t)). 


(6.  18) 


(6.  19) 


The  first  integral  in  (6.  18)  is  evidently  the  power  of  the  tractions 

external  with  respect  to  the  material  occupying  f)  in  the  reference  configu- 

+ 

ration;  the  second  integral,  which  vanishes  when  F(t)=F(t),  represents  the 
contribution  to  U'(t)  arising  from  the  deformation-gradient  discontinuity  at 
£{t).  Moreover,  P is  the  energy-momentum  tensor  originally  introduced 
by  Eshelby  [10]^  into  the  theory  of  defects  in  elastic  bodies.  Eshelby  es- 
tablished the  relation  between  the  tensor  P and  the  "force  exerted  on  the 
defect";  this  relation  — in  the  absence  of  such  defects  — gives  rise  to  a 
familiar  conservation  law  in  finite  and  linearized  ela stostatic s.  The  two- 
dimensional  ve rsion  of  the  conservation  law  to  which  we  are  alluding  asserts 

This  recent  expository  paper  contains  references  to  Eshelby' s related 
earlier  work. 
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that 


r s^P(F)N4i  ^0 

r 


(6.  20) 


for  all  plane  equilibrium  deformations  (of  a homogeneous  hyperelastic  sol- 
id) with  twice  continuously  differentiable  displacements,  if  T is  the  bound- 
ary of  an  arbitrary  regular  region  within  the  domain  of  regularity  of  the 
elastostatic  field  at  hand  and  N is  the  outward  unit  normal  vector  of  T.  The 
fracture-mechanical  implications  of  (6.  20)  and,  in  particular,  the  connec- 
tion between  the  so-called  ^-integral  and  the  energy  release  rate  at  the  tip 
of  a crack,  were  first  recognized  by  Rice  [11]. 

For  our  purposes  it  is  useful  to  cast  the  energy  identity  (6.  18)  into 
a different  form.  To  this  end  we  use  (6.  19),  (6.  7),  and  the  first  of  (6.  6)  to 
infer  that 


[p(t)j;M(t)  .L(t)  - -^(t)  .[F(t)];L(t)  -0, 

whence 

[P(t)]'N(t)  =[[P(t)]'N(t)  '^(t)}^ (t). 


Accordingly,  the  energy  identity  (6.  18)  is  equivalent  to 


where 


,(t) 


U’(t)=  £($- t)*v(^>t)c^ft-[H(t)]"JV(5,  t)  .N(t)d§, 


?j(t) 


H(t)-H(F(t),N(t))  , H(t)-H(F(t),N(t)). 


In  view  of  (6.  24)  and  (6.  7)  one  has 


(6.21) 


(6.  22) 


(6.23) 


M6.24) 


[H(t)]- =s(t)  -[F(t)];N(t)-[W(t)];  , 


(6.  25) 
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so  that  the  jump  in  H(t)  is  the  excess,  over  the  jump  in  the  strain- ene rgy 
density  W(t),  of  the  "internal  work"  done  by  the  Piola  traction  ^(t)  through 
the  jump  [F(t)]^^(t)  in  the  deformation  gradient  normal  to  the  material 
shock-line  £(t). 

Bearing  in  mind  the  physical  significaince  of  the  second  term  in 
(6.23)  as  the  contribution  to  the  energy  rate  U"  (t)  due  to  the  deformation- 
gradient  discontinuity  at  the  shock-line  £(t),  we  introduce  the  following  def- 
inition. ^ time -dependent  elastostatic  shock  is  dissipative  at  an  instant  t 
in  the  interior  of  its  interval  of  duration  [tj  ,^2^—  and  only  jd 

?2(t) 

[H(t)];Jv(§.t).N(t)d§>0  (6.26) 

?l(t) 

for  every  test  domain  £ admissible  at  this  instant. 

We  now  establish  necessary  and  sufficient  conditions  in  order  that 
a time-dependent  piecewise  homogeneous  elastostatic  shock  be  dissipative 
at  a given  instant.  Suppose  such  a shock  is  dissipative  at  t.  Then,  we 
show  first,  its  material  shock-line  £ must  be  in  a state  of  translation  at 
that  moment  and  the  direction  of  this  instantaneous  translation  cannot  be 
parallel  to  £(t): 

y(?,t)=V(0,t)HV(t)  (-oo<?<oo)  , V(t).N(t)^0.  (6.27) 

Indeed,  if  the  first  of  (6.  27)  were  false,  there  would  exist  an  instantaneous 
center  of  rotation  for  the  rigid  motion  of  £ at  the  time  t.  In  this  event 
V ( S,  ,t)*N(t)  would  range  continuously  over  all  real  numbers  as  § ranges 
over  (-03,00)  and  thus  the  inequality  (6.26)  could  not  possibly  hold  for  eve ry 
jft  tliat  is  the  interior  of  a sufficiently  small  circle  centered  on  £(t).  Further, 
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t 


the  first  of  (6.  27)  is  incompatible  with  (6.  26)  unless  the  second  of  (6.  27) 
holds  true.  Next,  from  (6.  26),  (6.  27)  follows 

[H(t)J'>0  (6.2«) 

provided  £(t)  is  oriented  by  means  of  L(t)  in  such  a way  that 

V(t).N(t)>0.  (6.29) 

Consequently,  H(t)  must  increase  as  £(t)  is  traversed  from  A(t)  to  ft(t),  if 
^(t)  is  understood  to  be  the  half-plane  into  which  Z is  advancing  at  the  in- 
stant under  consideration.  Also.  (6.  28),  (6.  29)  and  the  first  of  (6.  27)  are 
evidently  sufficient  for  the  dissipativity  of  the  shock  at  tinie  t.  Finally,  it 
is  clear  from  the  above  that  a piecewise  homogeneous  time -dependent  shock 
cannot  be  dissipative  throughout  a time-interval  unless  the  motion  of 

its  material  shock-line  is  purely  translatory. 

The  dissipation  inequality  (6.  28),  as  will  become  apparent  later  on, 
assumes  a role  analogous  to  that  played  by  the  entropy  inequality  in  connec- 
tion with  gas -dynamical  shocks.  Our  immediate  goal  is  to  deduce  a useful 
alternative  representation  for  the  jump  [H(t)]^  appropriate  to  a time-depend- 
ent elastostatic  shock.  The  subsequent  considerations  apply  to  any  fixed  in- 
stant in  the  time  interval  (t^.t^);  for  the  sake  of  brevity  we  shall  suppress 
the  argument  t in  the  equations  to  follow.  Because  of  (6.  24)  and  (3.  3),  (3.  32) 
we  are  entitled  to  write 

H(F,N)- -W(F)+J(F)FN.^(F)(F‘S^N  if  F=^  or  F-F,  (6.  30) 

where  jt(F)  is  the  actual  stress  field  associated  with  the  nominal  stress 
field  ja(F).  On  tlu'  other  hand,  (3.  11),  (3.  12),  (3.  13)  and  the  traction  con- 
tinuity condition  in  the  form  (3.  33)  easily  lead  from  (6.  30)  to 


1 
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H(F,N)=  -W(F)  + 


n "GIFjn 


+ 


(6.  31) 


in  which  n is  the  unit  normal  vector  of  the  (moving)  spatial  shock-line 
and^(F)^(F)  the  actual  traction  vector  along  Recall  now  that  £ is  the 

unit  direction  vector  of  and  X*  the  (time -dependent)  coordinate  frame 
with  the  base  vectors  (f  ,n)  and  the  same  origin  as  X.  On  expanding  the 
two  scalar  products  entering  (6.  31)  in  x‘ , noting  that  n and  t have  the  com- 
ponents (0,  1)  and  this  frame,  we  obtain 


/ / 1 ^ - 

H(g,X^)=  -WCg)+J(j^)  ^^-^(F)  + 

Finally,  from  (6.  32)  and  (3.  25),  (3.  27),  the  continuity  across  of  the 
stress  components  "^'12' ^Z2’  last  two  of  (6.24),  follows  the  useful 

re  suit 


[H];= 


[W]‘+J(HT'j2-eT22)  . e=l-6. 


(6.  33) 


where 


+/ 

■^22 


=^22(F)- 


(6.  34) 


In  (6.  33),  K and  6 are  the  values  of  the  amount  of  shear  and  of  the  relative 
stretch  at  right  angles  to  inherent  in  the  instantaneous  supplementary 
deformation.  ^ 

At  this  point  we  employ  (6.  33)  in  applying  the  dissipation  condition 
(6,28)  to  the  special  case  of  a piecewise  homogeneous  time -dependent  elas- 
tostatic  shock  occurring  in  Material  1.  To  this  end  we  observe  first  that 
the  results  of  Section  5 continue  to  hold  at  present  if  proper  allowance  is 
made  for  the  time-dependence  of  the  equilibrium  shock.  In  particular,  the 

^Refer  to  the  decomposition  (6.8)  and  to  (3. 19),  (3.28).  See  also  Figure  6. 
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jump  [VV(t)j"  in  tile  energy  density  is  available  from  (5.  22)  if  in  this  equa- 
tion and  6 are  replaced  by  \ ^ (t) , (f) . P(t)  ,cp(t) , and  6(t).  Fur- 


ther, (1.  14),  (1.  18)  anti  (4.41)  yield  the  stress  components 

,=U  -jl— on 


12'’ j .?^12 


22 


(6.  35) 


Tile  required  jump  in  H(t)  may  now  be  calculated  from  (6.  33)  with  the  aid  of 


+ / 


(5.22),  (6.  35),  the  formulas  (5.  4)  for  and  equation  (5.8)  for  the  amount 


of  supplementary  shear  k.  This  computation  ultimately  gives 

,2._  , 3 


[Hj;= 

2j(l-e)  (l+2pcos2cp+(3  ) 


(6.  36) 


Now  p >0 , J>0 , who  rea s ] P | < 1 and  e <1  acco rding  to  (5.6).  Hence  (6.  36 ) 
impl  ie  s 

(6.  37) 


fH(t)]^>t)  if  and  only  if  e(t)>0, 


so  that  the  dissipation  condition  (6.  28),  which  rests  on  the  assumption  that 

■(- 


the  shock  is  advancing  into  J<(t)  at  the  time  t,  leads  to  the  shock  - strength 
restriction  0(t)>O  in  the  case  of  Material  1 (for  all  admissible  homogeneous 
deformations  on  !ft(t)).  This  restriction  was  anticipated  in  Section  5,  when 
the  detailed  discussion  of  some  of  the  results  deduced  there  was  confined  to 
non-negative  values  of  ?. 

One  gathers  from  (6.  36)  that  for  Material  1 , 


ase^O 


(6. 38) 


at  any  particular  instant.  The  gi.^neral  validity  of  this  estimate,  even  for 
elastostatic  shocks  in  anisotropic  hyp^  relastn.  matiTials,  can  be  established 


J 
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by  means  of  an  appropriate  weak-shock  expansion^.  Furthermore,  (6.  38) 
is  the  analogue  of  a familiar  property  of  the  entropy  jump  in  gas -dynamical 
shocks . 

We  now  list  various  implications  of  the  conclusion  that  e(t)>0  for  a 
dissipative  shock  in  Material  1 if  the  shock  happens  to  be  advancing  into 
R(t)  at  this  moment.  The  following  results  apply  to  the  instant  under  con- 
sideration, although  the  argument  t will  be  suppressed. 

First,  in  view  of  the  kinematic  significance  of  the  shock- strength 
C,  the  supplementary  deformation  involves  a contraction  at  right  angles  to 
the  spatial  shock-line  Second,  it  is  clear  from  (5.  3)  and  (5. 6),  (5. 8) 

that  the  sign  of  the  supplementary  amount  of  shear  h depends  on  the  nature 

of  the  instantaneous  principal  stretches  inherent  in  the  given  homogeneous 
+ 

deformation  on  R . One  has 

M.<0  if  (3>0  (X^>Xj)  , K>0  if  (3<0  (:^2<XJ)  (6.39) 

for  shocks  with  a spatial  shock-angle  in  the  first  quadrant  (0,tt/2). 

2 

Next,  from  (4.4)  , (5.6),  and  (5.22)  one  draw's 

[J]+<0  . [P]>0  , [W]'<0.  (6.40) 

Consequently,  in  the  present  instance,  the  area-ratio  and  the  energy 
density  decrease,  while  the  mass  density  increases  — as  the  material 
shock-line  £ is  traversed  from  R to  R.  Further,  we  recall  from  the  discus- 
sion at  the  end  of  Section  5 that  in  the  case  of  Material  1,  e>0  assures  a 

4 

loss  of  ellipticity  on  R at  all  admissible  deformations  on  in  contrast, 

^ k'or  this  purpose  one  requires  a weak -shock  expansion  of  higher  order 
than  that  considered  in  Section  4. 

2 

Note  that  (4.4)  holds  true  for  shocks  of  finite  strength. 


the  ensuing  deformation  on  R may  or  may  not  be  accompanied  by  a failure 
of  ellipticity  of  the  displacement  equations  of  equilibrium.  Thus,  at  least 
for  Material  1,  the  situation  is  parallel  to  that  encountered  in  gas -dynamical 
shocks,  where  — as  a consequence  of  the  entropy  inequality  — the  flovv"  is 
necessarily  supersonic  on  the  upstream  side  but  may  be  either  supersonic 
or  subsonic  on  the  downstream  side. 

Finally,  consider  a time -dependent  piecewise  homogeneous  equilib- 
rium shock  in  Material  1 that  is  dissipative  at  all  times  in  the  interior  of 
its  interval  of  duration  [t^,t2].  Suppose,  in  addition,  that  the  material 
shock-line  Z(t)  of  the  shock  advances  steadily  into  ft(t)  for  t^<t<t2-  Then, 
e>0  on  (tj^,t2)  and  — as  was  shown  earlier  in  this  section  — the  motion  of 
£ is  one  of  pure  translation  throughout  this  range  of  time.  Thus,  the 
material  shock-angle  obeys 

<J>(t)=$o  (6.41) 

in  which  is  a constant.  In  these  circumstances  according  to  (3.  14), 
for  a given  (admissible)  deformation  on  R(t)  (t^st^t^),  the  spatial  shock- 
angle  cp  becomes  a known  function  of  <i>^  and  the  time,  i.  e.  , 

Cp=Cp(t;$^)  confirms  easily  that  the  first  of  (5.9)  may 

now  be  inverted  to  give  e=e(t;4>^),  so  that  the  shock -strength  becomes  a 
fully  determinate  function  of  time  for  every  fixed  value  of  We  shall 

not  explore  the  specific  properties  of  e(t;«I>^),  however,  since  the  fore- 
going considerations  arc  strictly  liniited  to  piecewise  homogeneous  shocks 
and  have  no  local  analogue  in  connection  with  more  general  equilibrium 
shocks . 

For  comparison  purposes  we  cite  also,  without  proof,  some  results 
obtained  by  applying  the  dissipation  condition  (6. 28)  to  a weak 
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time -dependent  shock  in  Material  2.  Here  one  finds  that  e<0,  so  that  the 
instantaneous  supplementary  deformation  involves  a stretching  perpen- 
dicular to  Z^.  On  the  other  hand,  (4.42)  and  (4.  36)  imply  that  the  in- 
equalities (6.  39),  governing  the  sign  of  the  supplementary  amount  of  shear 

K,  remain  valid  for  Material  2 if  |e]  is  sufficiently  small.  Also,  since  j 

i 

e<0,  the  first  two  inequalities  in  (6.40)  now  give  way  to  j 

t 

[J]‘>0  , [p]‘<0.  (6.42) 

The  sign  of  the  energy-jump  [W]^  can  no  longer  be  inferred  from  the  sign 
of  e alone  in  this  special  case;  it  is  found  to  depend  on  the  instantaneous 
character  of  the  deformation  pre-assigned  on  R. 

7 . Equilibrium -shock  formation  as  a bifurcation  process . Luders  bands. 

Discussion. 

The  analysis  in  Section  4 led  to  the  conclusion  — valid  under  certain 
assumptions  spelled  out  there  — that  the  existence  of  a piecewise  homo- 
geneous elastostatic  shock,  and  hence  the  emergence  of  discontinuities  in 
the  first  deformation  gradients,  is  contingent  upon  a breakdown  of  ellip- 
ticity  in  the  displacement  equations  of  equilibrium  associated  with  the 
homogeneous  deformation  prescribed  on  R.  As  is  well  known  and  easily 
verifiable,  this  failure  of  ellipticity  in  turn  renders  the  given  deformation 
dynamically  unstable:  there  are  initially  periodic,  small-amplitude,  dis- 

turbances of  the  above  uniform  elastostatic  field  that  give  rise  to  solutions 
of  the  appropriate  linearized  displacement  equations  of  motion  which  grow 
beyond  bounds  with  time . 


The  foregoing  state  of  affairs  suggests  the  possibility  of  viewing 
the  orocess  of  shock  formation  as  a bifurcation  from  a homogeneous 
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equilibrium  field  in  an  elastic  body.  Such  an  interpretation  of  equilibrium 
shocks  leads  one  to  wonder  about  their  relevance  to  the  familiar  failure 
phenomenon  of  Luders  bands,  commonly  observed  in  mild  steels,  which 

has  received  repeated  analytical  attention  within  plasticity  theory.  ^ 

..  2 

Luders  bands  are  known  to  develop  in  specimens  subjected  to  uni- 
axial tension  or  compression  once  the  loads  have  exceeded  the  yield  limit; 
their  inclination  relative  to  the  load  axis  is  often  close  to  45  degrees,  al- 
though considerable  departures  from  this  angle  have  been  reported.  Fig- 
ure 10,  which  follows  Nadai  [13],  shows  a schematic  diagram  of  Luders 
bands.  The  experimental  findings  indicate  abrupt  changes  in  the  deforma- 
tion giadients  across  the  interfaces  between  such  a band  and  the  adjacent 
material,  the  deformation  within  each  band  being  predominantly  one  of 
shear  parallel  to  the  interfaces  Further,  the  sense  of  this  shear  deforma- 
tion undergoes  a reversal  as  the  loading  is  changed  from  tension  to  com- 
pression. 

In  assessing  the  extent  to  which  predictions  based  on  the  theory  of 
piecewise  homogeneous  equilibrium  shocks  resemble  experimental  observa- 
tions pertaining  to  the  formation  of  Luders  bands,  it  should  be  kept  in  mind 
that  the  actual  test  situation,  which  is  in  fact  three-dimensional,  comes 
closer  to  conditions  of  plane  stress  than  to  plane  strain.  Moreover, 
elementary  equilibrium  considerations  preclude  the  existence  of  a homo- 
geneous field  of  deformation  and  stress  in  a Luders  band  of  finite  width 
(see  Figure  10)  since  such  a uniform  field  is  incompatible  with  the  boundary 

^See,  for  example,  Thomas  [12] 

2 

See  Nadai  [13],  Chapter  18,  for  a detailed  description  of  the  pertinent 
experimental  observations  and  for  references  to  the  previous  phenomeno- 
logical literature  on  Luders  bands 
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conditions  at  the  load-free  parallel  edges  of  the  specimen.  For  this  reason 
a piecewise  homogeneous  elastostatic  shock  cannot  possibly  describe  the 
observed  behavior  near  the  edges  of  the  slab.  On  the  other  hand,  it  is  a 
simple  matter  to  symmetrize  the  elastostatic  field  of  an  equilibrium  shock 
and  thus  arrive  at  a piecewise  homogeneous  field  involving  three  homo- 
geneous zones:  two  half-planes  in  which  the  same  given  deformation  is 
sustained,  separated  by  a strip  that  undergoes  a distinct  deformation,  the 
displacements  and  cractions  being  continuous  across  the  interfaces. 

Despite  the  limitations  pointed  out  above,  the  present  theory  of 
equilibrium  shocks  exhibits  certain  striking  features  that  support  its 
relevance  (within  the  framework  of  continuum  mechanics)  as  far  as  Lliders 
bands  are  concerned.  The  following  considerations  pertain  to  the  local 
behavior  near  a point  such  as  P in  Figure  10,  on  an  interface  of  a Lliders 
band. 

To  fix  ideas  and  solely  for  illustrative  purposes  we  draw  once 
more  on  the  special  isotropic  hyperelastic  Material  1.  Suppose  an  all- 
around  infinite  slab  of  this  material  occupies  the  entire  (x^  , x^) -plane  R 
in  its  undeformed  configuration  and  is  subjected  to  the  time -dependent 
quasi -static  pure  homogeneous  plane  deformation: 

j^=^(x,  t)  = A(t)x  for  all  x in  R (0st<oo),  (7.  1) 

where 

^^(t)  0 ■ 

A(t)=  (0st<oo),  (7.2) 

. 0 \2(t)_ 


A(t)  being  the  deformation -gradient  tensor  and  \Jt)  the  principal  stretches 
at  time  t.  We  shall  assume  that  the  foregoing  deformation  process  starts 


from  the  undeformed  state  at  the  instant  t=0  and  at  all  times  thereafter 
corresponds  to  a state  of  uni-axial  tension  or  compression  parallel  to  the 
X2*a-xis.  Thus}  x 


XJ0)  = 1 . T^^(t)=o^^(t)=0  . \j(t)=X'"(t), 
T22(t)=p[l-X2'®^®(t)]  , a22(t)=p[X2’^(t)-X2^  (t)]  (Ost<oo). 


(7.  3) 


In  addition  we  require  that  in  the  case  of  uni-axial  tension. 


X2(t)  is  steadily  increasing,  X2(t)"*oo  as  t-»oo. 


while  in  the  case  uni-axial  compression 


X2(t)  is  steadily  decreasing,  X2(i)"^^  t-^oo. 


(7.5) 


We  now  recall  from  Section  2 (see  Figure  1)  that  the  homogeneous 
time -dependent  equilibrium  deformation  characterized  by  (7.  l),  (7.2), 

(7.  3)  entails  a loss  of  ellipticity  whenever  the  stretch  X2(t)  equals  or 
exceeds  a certain  critical  value  in  the  tension  case  or  fails  to  be  above 
another  critical  value  in  the  compression  case.  Hence,  by  virtue  of  (7.4), 
(7.  5),  in  either  case  there  is  an  instant  t^  such  that  the  deformation  pro- 
cess under  consideration  is  accompanied  by  a loss  of  ellipticity  — and 
thus  also  by  a loss  of  dynamic  stability  — at  all  times  t^t^.  Moreover,  it 

is  now  clear  from  the  conclusions  reached  in  Sections  5,  6 that  if  t.st 

1 o 

the  following  bifurcation  of  the  homogeneous  time -dependent  deformation 
at  hand  into  a piecewise  homogeneous  time -dependent  equilibrium  shock 
becomes  possible: 

^Recall  (2.8),  which  apply  to  uni-axial  stress  parallel  to  the  Xj-axis. 
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provided 


F(t)x  for  all  X inR(t),  F(l)=A(t) 
F(t)x+b(t)  for  all  x in  rt(t)  (t^st<'cx>), 

V(t)*N(t)>0  (t  st<co). 

~ ^ '1 


(7.6) 


:7.  7) 


Here  V(t)  is  the  velocity  of  the  material  shock-line  £(t)  dividing  into  the  two 
+ " 

half-planes  R(t),  R(t);  also  £(t),  whose  unit  normal  vector  N(t)  points  into 
ft(t),  is  in  translation  because  of  the  dissipativity  of  the  shock.  The  assump- 
tion (7.  7)  amounts  to  the  physically  motivated  requirement  that  the  deforma- 
tion in  R(t)  — associated  with  the  local  deformation  in  the  Lliders  band  — stead- 
ily encroaches  upon  the  as  yet  unencumbered  deformation  prevailing  on  ^(t). 
Because  of  (7.  7)  we  may  appeal  to  (6. 37)  and  the  dissipativity  of  the  shock 
to  infer  that  its  strength  e(t)>0  for  t^<t<oo. 

The  results  pertaining  to  Material  1 in  Sections  2,4,6  now  permit 
various  inferences  that  bear  on  the  significance  of  equilibrium  shocks  in 
connection  with  Luders  bands.  First,  in  the  tension  case  the  shock- 
bifurcation  (7.6)  can  arise  only  after  the  nominal  stress  022  associated 
with  the  homogeneous  deformation  (7.  1)  has  passed  its  peak  (see  Figure  1). 
Second,  in  both  the  tension  and  the  compression  case,  the  angle  of  inclina- 
tion of  the  ''emerging''  spatial  shock-line  relative  to  the  load-axis,  i.e. 
the  limiting  value  of  this  angle  at  zero  shock-strength  and  at  an  incipient 
failure  of  ellipticity,  is  45  degrees  for  Material  1.  In  general  this  angle 
is  a material  property  and  may  be  different  for  tension  and  compression. 

Next,  the  shock -formation  in  the  present  instance  involves  a sup- 
plementary contraction  at  right  angk's  to  the  interface  between  the  half- 
planes  separated  by  the  spatial  shock-line.  On  the  other  hand,  according 


A 


to  (6.  59),  the  amount  of  shear  h accompanying  each  instantaneous  supple- 
mentary deformation  is  positive  ^ in  the  case  of  a con^pression -induced 
shock  and  negative  for  a tension-induced  shock  in  Material  1.  The  last- 
mentioned  conclusion  is  at  variance  with  the  usually  observed  deformation 
pattern  of  Lilders  bands,  schematically  depicted  in  Figure  10.  We  do  not 
know  whether  the  opposite  sense  of  over-all  shear  has  been  encountered  in 
tests  of  actual  materials.  Be  this  as  it  may,  one  can  show  by  asymptotic 
means  that  there  are  strain-ent^rgy  densities,  even  within  the  limited 
class  characterized  by  (2.  1),  that  lead  to  h>0  in  weak  tensile  - shocks . 

The  above  predictions  concerning  the  sign  of  k therefore  merely  reflect  a 

2 

peculiarity  of  the  particular  idealized  material  under  discussion  . 

The  preceding  remarks  suggest  a word  of  caution  regarding  the 
highly  special  nature  of  the  two  hypothetical  materials  used  in  this  paper 
to  illustrate  the  theory  of  elastostatic  shocks.  With  a view  toward  physical 
ly  realistic  applications  of  the  theory  it  would  seem  essential  to  explore  a 
v/ider  range  of  liype  relasti c solids  that  can  sustain  a loss  of  ellipticity. 

In  particular,  it  would  be  of  interest  to  construct  ideal  materials  that 
adtiiit  a loss  of  ellipticity  in  uni -axial  tension  and  compression  at  pre- 
assignable  stretch  levcTs  Unfortunately,  such  an  adaptability  of  the  as- 
sumerl  constitutive"  behavior  is  bound  to  inc\ir  mounting  analytical 
complexities. 

Tht’  analysis  carried  o\i1  in  this  investigation  may  Ije  generalized 

^See  Figure  6(b)  for  ti^'-  geomi-tric  meaning  of  >'.>0. 

weak  tensile  shoik  in  Material  2 also  involvi-s  i<0,  as  follows  from  the 
observatif)ns  at  the  end  of  Section  tj . We  recall  that  a compressive  shock 
cannot  occur  in  Material  2 sinci'  the  latter  dfies  not  suffer  a loss  of  ellipti- 
city in  uni -axial  compression. 
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to  encompass  no  longer  piecewise  homogeneous  plane  elastostatic  deforma- 
tions with  continuous  displacement  fields  and  merely  piecewise  continuous 
first  and  second  displacement  gradients,  whose  discontinuities  are  permit- 
ted to  occur  along  curved  shock-lines.  Further,  it  is  a priori  clear  that 
the  results  deduced  here  for  piecewise  homogeneous  equilibrium  shocks 
at  once  apply  to  the  local  situation  at  an  interior  point  of  a curved  shock- 
line in  a non-homogeneous  elastostatic  field.  The  emergence  of  weak 
shocks  of  this  more  general  type  evidently  necessitates  a breakdown  of 
ellipticity  in  the  displacement  equations  of  equilibrium  at  the  shock-line, 
which  — in  the  weak-shock  liniit  — must  be  a characteristic  line  associated 
with  these  equations  at  the  prevailing  deformation.  Also,  the  dissipation 
inequality  (6.28)  remains  locally  valid  in  the  present  circumstances,  but 
the  dissipativity  of  such  a shock  no  longer  requires  that  the  motion  of  the 
shock-line  be  translatory. 

Presumably  non-homogeneous  equilibrium  shocks  would  arise  in 
boundary- value  problems  of  finite  elastostatics  at  loads  that  cause  suf- 
ficiently severe  local  deformations  — severe  enough  to  induce  a local 
failure  of  ellipticity  in  the  governing  field  equations,  provided  of  course 
the  underlying  elastic  potential  admits  such  a failure^.  In  the  presence 
of  geometric  or  material  sources  of  stress  concentrations  (such  as  holes, 
notches,  or  inclusions)  elastostatic  shocks  could  thus  evolve  at  compara- 
tively moderate  loads. 


^ The  situation  envisaged  here  is  analogous  to  that  encountered  in  boundary- 
value  problems  for  stationary  flows  in  gas  dynamics  that  give  rise  to  super- 
sonic regions  and  concomitant  shocks. 
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FIGURE  I.  MATERIAL  I.  PLANE-STRAIN  UNI-AXIAL  STRESS:  T2  = 0. 
NOMINAL  AND  ACTUAL  STRESS  VS  AXIAL  STRETCH 
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FIGURE  3.  MATERIAL  2.  PLANE-STRAIN  UNI-AXIAL  STRESS:  T2  =0. 
NOMINAL  AND  ACTUAL  STRESS  VS  AXIAL  STRETCH 
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FIGURE  5.  KINEMATICS  OF  PIECEWISE  HOMOGENEOUS  ELASTOSTATIC  SHOCKS 
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This  investigation  concerns  equilibrium  fields  with  discontinuous  displace- 
ment gradients,  but  continuous  displacements,  in  the  theory  of  finite  plane  deforma- 
tions of  possibly  anisotropic,  compressible  elastic  solids.  d'Elastostatic  shocks""  of 
this  kind,  which  resemble  in  many  respects  gas -dynamical  shocks  associated  with 
steady  flov/s,  are  shown  to  exist  only  if  and  when  the  governing  field  equations  of 
equilibrium  suffer  a loss  of  ellipticity.  The  local  structure  of  such  shocks,  near 
a point  on  the  shock -line,  is  studied  with  particular  attention  to  weak  shocks,  and 
an  example  pertaining  to  a shock  of  finite  strength  is  explored  in  detail.  Also,  neces 
sary  and  sufficient  conditions  for  the^"dissipativity""  of  time -dependent  equilibrium 
shocks  are  established.  Finally,  the  relevance  of  the  analysis  carried  out  here  to 
localized  shear  failures  — such  as  those  involved  in  the  formation  of  Liiders  bands  — 
is  discussed. 
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